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Abstract 

An  experimental  investigation  of  two-dimensional 
confined  Jet  thrust  vector  control  nozzles  was  performed. 
Thrust  vector  control  was  accomplished  by  using  secondary 
flow  injection  in  the  divergent  section  of  the  nozzle. 
Schlieren  photographs  and  video  tapes  were  used  to  study 
flow  separation  and  internal  shock  structures.  Nozzle 
performance  parameters  were  determined  for  nozzle  flow  with 
and  without  secondary  flows.  These  parameters  included 
nozzle  forces,  vector  angles,  thrust  efficiencies,  and  flow 
switching  response  times.  Vector  angles  as  great  as  18 
degrees  with  thrust  efficiencies  of  0.79  were  measured. 
Several  confined  Jet  nozzles  with  variations  in  secondary 
flow  port  designs  were  tested  and  results  were  compared  to 
each  other.  Converging -diverging  nozzles  of  similar  design 
to  the  confined  Jet  nozzles  were  also  tested  and  results 
were  compared  to  the  confined  Jet  nozzle  results. 

Existing  prediction  models  for  nozzle  side  to  axial 
force  ratio  were  evaluated.  A  model  for  nozzle  total  forces 
based  on  shock  losses  that  predicted  values  very  close  to 
actual  results  was  developed. 
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TWO-DIMENSIONAL  CONFINED  JET  THRUST  VECTOR  CONTROL: 


OPERATING  MECHANISMS  AND  PERFORMANCE 


I .  I ntroducti on 


Background 

An  important  part  of  any  airborne  or  space  vehicle  is 
its  attitude  and  guidance  control  system.  Most  aircraft  and 
missiles  use  aerodynamic  surfaces  to  adjust  their  attitude 
and  direction  of  flight.  This  type  of  control  varies 
greatly  with  the  speed  and  altitude  of  the  vehicle  and  has 
the  disadvantage  of  added  drag.  Changing  the  direction  of  a 
vehicle’s  thrust  is  a  method  of  augmenting  or  ever,  replacing 
aerodynamic  control  surfaces. 

Possible  methods  of  thrust  vector  control  C  TVO  for 
rocket  engines  can  be  divided  into  three  groups.  Group  1 
are  those  nozzles  which  use  direct  mechanical  interference 
with  the  exhaust  gases  such  as  jet  vanes,  spoilers,  and 
jetavators.  Group  2  systems  require  the  movement  of  the 
whole  thrust-producing  apparatus  such  as  swivel  or  gimballed 
nozzles.  Both  of  these  groups  require  high-temperature 
materials  and  elaborate  seals.  Group  3  systems  use 
nonmechanical  interference  with  exhaust  gases  such  as  gas 
bleeding  or  the  secondary  injection  of  gas  into  the  nozzle 
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C  4 :  2D  .  Gas  bleeding  is  simply  the  release  of  combustion 
gases  through  a  channel  at  some  angle  to  the  centerline  axis 
of  the  nozzle.  Its  effective  side  force  is  therefore 
limited  to  the  momentum  of  the  gas  bleed  flow.  The 
secondary  injection  of  gas  into  the  nozzle  produces  a  side 
force  that  is  a  combination  of  the  secondary  jet  momentum 
and  the  interaction  Cdynamic  and  chemicalD  of  the  secondary 
jet  with  the  primary  axial  flow  C2:214D.  This  interaction 
allows  the  side  force  to  be  greater  than  the  momentum  of  the 
secondary  jet. 

Secondary  In  lection  Thrust  Vector  Control  CSITVCD 

Boundary  Layer  Thrust  Vector  Control  C BLTVCD  is  one  type 
of  SITVC.  It  uses  an  over -expanded ,  conver gi ng-di ver gi ng 
nozzle  with  control  ports  in  the  diverging  section  to  allow 
ambient  air  as  a  secondary  gas  into  the  nozzle.  When  a  port 
is  opened,  the  pressure  rise  due  to  the  ambient  air  augment  ; 
the  exit ing  separation  region  in  the  nozzle.  This  causes 
the  primary  flow  path  to  attach  to  the  opposite  wall  where 
the  pressure  is  lower.  Figure  1  illustrates  this  flow 
mechanism  This  type  of  TVC  requires  that  the  ambient 
pressure  be  greater  than  the  pressure  at  the  secondary  port 
in  the  nozzle.  BLTVC  systems  will  therefore  have  altitude 
limitations  to  their  operational  envelope  C  8:1, 2D. 

Confined  Jet  Thrust  Vector  Control  CCJTVCD  is  a  type  of 
ST  TVC  that  may  not  have  the  altitude  restrictions  of  BLTVC. 
It  uses  a  reconverging  section  on  a  nozzle  to  contain  the 
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Fig  1  Boundary  Layer  Thrust  Vector  Control 
CBLTVO  System 
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Fig  2  Confined  Jet  Thrust  Vector  Control  CCJTVC)  System 


separated  region,  thereby  insulating  it  from  ambient 
conditions.  Secondary  flow  enters  the  nozzle  from  some 
source  and  deflects  the  axial  flow  to  the  opposite  wall, 
where  it  attaches.  Figure  2  illustrates  a  typical  CJTVC 
system.  Besides  its  increased  operational  envelope,  CJTVC 
usually  has  relatively  small  secondary  mass  flow  rates. 

This  means  that  the  secondary  flow  supply  storage 
requirement  is  relatively  small  and  can  be  readily  packaged 
for  typical  rocket  TVC  applications.  In  addition  to  the 
compressed  secondary  gas  supply,  small,  fast-responding 
control  valves  are  required  for  the  system  C8:8J. 

Previous  CJTVC  Work 

Several  CJTVC  designs  have  been  constructed  and  tested 
at  the  Air  Force  Institute  of  Technology  CAFITJ.  Porzio 
designed  and  tested  an  axl symmetric  CJTVC  nozzle  0  2:7} 
variations  of  which  other  students  tested.  Cates  designed 
and  tested  several  two-dimensional  CJTVC  nozzle 
configurations,  only  one  of  which  could  be  consistently 
vectored  at  his  test  conditions  C5:47}.  Talda  improved 
Cates  design  and  developed  several  CJTVC  configurations 
which  could  consistently  vector.  Talda’s  nozzle  has 
efficiencies  as  high  as  855-#,  with  ±  20  degrees  deflections 
with  vectoring  response  times  as  low  as  50  ms  C16:104J. 
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The  goal  of  this  thesis  is  to  continue  the  previous  work 
in  the  design  and  testing  of  two-dimensional  CJTVC  nozzles. 
Emphasis  is  placed  on  studying  existing  designs  to  determine 
their  thrust  vector  mechanisms  in  hopes  of  developing  nozzle 
design  optimization  parameters.  Specific  objectives  for 
reaching  this  goal  are: 

1.  Use  flow  visualization  techniques  to  determine  flow 
separation  points  within  the  nozzles. 

2.  Use  flow  visualization  techniques  to  study  shock 
development  and  interaction  within  the  nozzles. 

3.  Measure  and  analyze  nozzle  performance 
characteristics  at  varying  operating  conditions.  These 
performance  characteristics  consist  of  side  and  axial  force 
measurements,  flow  vector  angles,  nozzle  thrust 
efficiencies,  and  side  momentum  gain.  The  momentum  gain  is 
the  ratio  of  the  measured  side  force  to  the  momentum  of  the 
secondary  flow.  Operating  conditions  will  be  varied  by 
changing  the  nozzle  pressure  ratio  CNPR3  and  the  secondary 
to  primary  mass  flow  ratio  CMFR3.  NPR  is  the  ratio  of 
primary  flow  total  pressure  to  ambient  pressure.  MFR  is  the 
ratio  of  the  secondary  mass  flow  rate  to  the  primary  mass 
flow  rate 

4.  Study  nozzle  vector  switching  characteristics, 
concentrating  on  determining  switching  response  time  and 
stable  operating  conditions.  Switching  response  time  is  the 
time  it  takes  for  primary  flow  to  vector  measured  from  the 
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point  when  secondary  flow  is  introduced. 


5.  Develop  prediction  models  for 
particularly,  a  model  that  can  predic 
occur  in  the  nozzle.  Evaluate  sever a 
prediction  models. 


nozzle  performance, 
the  losses  which 
existing  S1TVC 
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Talda's  LM  configuration  had  the  best  vectoring 
performance  of  all  the  two-dimensional  CJTVC  nozzles 
reviewed.  It  is  a  modification  to  Cates  best  performance 
design  where  sharp  linear  walls  are  replaced  by  smooth 
curved  walls  in  the  reconverging  section  C16:  1£-16, 32}  . 
Therefore,  it  was  used  as  the  basic  nozzle  design  for  this 
study.  Talda’s  LM  nozzle  was  tested  at  a  maximum  nozzle 
pressure  ratio  C  the  ratio  of  primary  flow  total  pressure  to 
ambient  static  pressure)  of  approximately  15.  Typical  solid 
rockets  have  nozzle  pressure  ratios  CNPR}  ranging  from 
50-1000  C7:40}.  Testing  a  nozzle  at  these  higher  NPRs 
requires  either  extremely  high  primary  pressures  or 
extremely  low  back  pressures.  For  safety  and  convenience, 
this  study  used  reduced  back  pressure  chambers  to  increase 
test  NPRs.  A  one-third  scale  Clinear  dimensions}  version  of 
Talda’s  design  was  used  as  the  basis  for  all  nozzles  tested 
Cl  6: 13}.  The  reduced  size  nozzle  was  chosen  in  order  to 
reduce  mass  flow  rate  and  thereby  decrease  the  evacuation 
requirements  of  the  reduced  back  pressure  chambers.  Figure 
3  illustrates  the  nozzle  design  and  contains  relevant  nozzle 
geometry  data. 

A  prototype  scaled-down  nozzle  was  constructed  and 
tested  in  a  blowdown  wind  tunnel.  These  tests  confirmed  the 
vectoring  ability  of  the  nozzle  and  provided  some  flow 
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visualization  information.  Four  variations  of  the  Talda 
design  were  manufactured  for  testing: 

Nozzle  0:  Smooth  wall  configurations  Cno  secondary 

ports) 

Nozzle  1:  Two  secondary  ports  with  flow  perpendicular 
to  primary  flow  axis;  ports  located  three  throat  widths  from 
nozzle  throat 

Nozzle  2:  Same  as  nozzle  1  except  secondary  flow  is 
perpendicular  to  nozzle  wall  at  secondary  port  location 

Nozzle  3:  Four  SI  ports  with  flow  perpendicular  to 
primary  flow  axis;  two  ports  located  three  throat  widths 
from  nozzle  throat;  two  ports  located  at  maximum  nozzle 
area. 

In  addition  to  these  four  CJTVC  nozzles,  two 
converging-diverging  configurations  were  tested. 

Nozzle  4:  Exit  to  throat  area  ratio  of  9.5  Cbased  on 
the  maximum  area  to  throat  area  ratio  of  the  Talda  nozzle 
design) 

Nozzle  5:  Exit  to  throat  area  ratio  of  4.0  Cbased  on 
the  exit  to  throat  area  ratio  of  the  Talda  nozzle  design} 

Figure  4  provides  a  comparison  among  all  six  test 
nozzl es . 

Test  Model  Design 

All  previous  CJTVC  studies  at  AFIT  tested  the  nozzles 
using  static  test  stands  with  the  nozzles  attached  to  a 
force  measurement  load  cell.  Since  installed  performance  of 
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a  TVC  nozzle  is  important  in  evaluating  its  practical 
applications,  a  wind  tunnel  model  capable  of  supporting  the 
test  nozzles  was  designed,  fabricated,  and  tested.  Figure  5 
illustrates  this  design.  The  design  point  of  the  model  was 
to  operate  in  a  Mach  2.2  wind  tunnel.  The  main  components 
of  the  model  are  listed  in  the  following  paragraphs. 

1.  Cylindrical  main  body  with  a  30°  half -angle  nose 
cone  to  properly  deflect  oblique  shocks  away  from  the  test 
area.  Body  is  hollow  with  enough  space  to  accommodate  a 
model  S  scani valve  C48  pressure  ports}. 

2.  Interchangeable  rear  section  of  main  body  that  can 
accommodate  either  a  nozzle  test  assembly  with  static 
pressure  measurement  ports  or  a  nozzle  test  assembly  for 
schlieren  photographs.  Because  its  cross  section  does  not 
match  the  cylindrical  main  body,  the  schlieren  assembly 
cannot  be  used  in  the  supersonic  free  stream. 

3.  Nozzle  test  assembly  with  36  static  pressure 
measurement  ports.  The  pressures  at  these  ports  were 
measured  using  a  scani valve  connected  to  an  ENDEVCO 
differential  pressure  transducer.  Figure  6  shows  the 
components  of  the  assembly,  including  the  distribution  of 
the  pressure  measurement  ports  across  the  nozzle. 

4.  A  rearward-swept  support  strut  that  interfaces  to  a 
six-component  force  balance.  Its  cross  section  is  a 
modified  diamond  airfoil  with  a  16*  half -angle  leading  and 
trailing  edge  to  properly  deflect  oblique  shocks  away  from 
the  test  area. 
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Appendix  A  contains  more  detailed  drawings  and 
dimensions  of  the  test  model. 

Test  Facilities 

The  primary  test  facility  for  this  study  was  the  S-l 
supersonic/transonic  wind  tunnel  located  at  the  von  Karman 
Institute  of  Fluid  Dynamics  CVKI5.  Rhode-St-Genese.  Belgium. 
The  S-l  wind  tunnel  is  a  continuous  closed  circuit  facility 
driven  by  a  615  kW  axial  flow  compressor.  A  40  cm  x  40  cm 
test  section  was  used  for  testing  where  supersonic  flows  of 
Mach  2.0  or  2.25  and  transonic  flows  between  Mach  0.4  and 
1.05  are  possible.  The  test  section  contains  a 
three-degr ee-of -freedom  traversing  mechanism  for  model 
support  as  well  as  variable  incident  mechanism.  The  tunnel 
is  also  equipped  with  a  schlieren  photograph  system  C18:233. 

Data  acquisition  was  performed  using  a  variety  of 
instruments.  A  summary  of  the  main  subsystems  used  in  data 
acquisition  follows: 

1.  Commodore  Extension  Interface  CCEI5  data  acquisition 
system  containing  a  12-bit  A/D  converter  with  16 
single-ended  input  channels  with  a  maximum  conversion  rate 
of  25  kHz.  The  CEI  interfaced  to  a  8032  PET-Commodor e 
computer  which  executed  the  appropriate  data  acquisition 
software. 

2.  VKI  6-component  force  balance  with  its  supporting 
amplifiers,  filters,  Wheatstone  bridges,  and  CEI  connection. 
Only  the  two  normal  force  channels  and  the  axial  force 


14 


channel  were  used.  Design  details  of  the  balance  are 
included  in  Appendix  A. 

3.  Model  48S3  scani valve  with  its  supporting  DC  motor 
drive,  power  supply,  position  indicator,  and  CEI  connection. 

4.  ENDEVCO  Model  851 0B-50  pi ezoresi sti ve  pressure 
transducer  with  its  supporting  ampl i f i er /f i 1  ter .  It  is  a 
differential  pressure  transducer  with  a  range  of  0-50  psig. 

5.  TEKTRONIX  Model  2430  digital  oscilloscope  capable  of 
simultaneous  acquisition  of  two  input  channels.  It  has  a 
maximum  digitizing  rate  of  100  megasamples  per  second  and  a 
real-time  useful  storage  bandwidth  of  40  MHz. 

6.  Hewlett-Packard  Model  3582A  spectrum  analyzer 
capable  of  single  or  dual  channel  sampling  and  digitizing. 

?.  Solenoid  control  panel  for  secondary  Injection 
valves  capable  of  sending  an  acquisition  trigger  signal 
simultaneous  to  the  valve  open/close  command  to  either  the 
digital  oscilloscope  or  spectrum  analyzer. 

8.  ROTAMETER  Model  G1  .  1000  mass  flow  meter.  Primary 
flow  rates  were  measured  using  a  steel  float  while  secondary 
flow  rates  were  measured  using  a  plastic  float. 

Prior  to  testing  at  VKI ,  the  wind  tunnel  model  was 
static  tested  at  facilities  in  building  454  of  Wright 
Patterson  AFB,  Ohio.  These  were  the  same  facilities  used  by 
Cates  and  Talda  for  their  CJTVC  nozzles.  The  tests 
confirmed  the  nozzles  vectoring  capability  and  structural 
integrity  prior  to  shipment  to  VKI.  Appendix  A  includes  a 
drawing  of  the  VKI  tunnel  and  how  the  instrumentation  was 
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integrated  into  it. 

The  test  model  met  the  weight  and  cross  sectional  area 
blockage  requirements  of  the  VKI  tunnel.  However,  when 
placed  on  the  tunnel  sting,  the  model  created  a  pitching 
moment  at  operating  conditions  that  would  damage  the  force 
balance.  Therefore,  force  measurements  were  not  possible 
with  the  wind  tunnel  model.  Using  a  solid  steel  ■'dummy" 
balance  Cno  strain  gages} ,  the  model  was  able  to  support  the 
schlieren  apparatus  for  flow  visualization  tests.  Force 
measurements  were  taken  by  mounting  the  nozzle  assembly 
directly  to  the  force  balance  as  previously  illustrated  in 
Fig  6.  Because  of  this  arrangement,  force  measurements 
taken  in  a  supersonic  freestream  were  not  possible. 
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Separation  Studies  Cno  secondary  injection!) 

CJTVC  nozzles  depend  on  the  primary  flow  separating  from 
one  wall  and  becoming  attached  to  the  opposite  wall  for  flow 
vectoring.  Therefore,  it  is  important  to  know  where  the 
flow  will  inherently  separate  in  the  nozzle.  The  distance 
along  the  nozzle  axis  to  where  the  flow  separates  from  the 
walls  for  nozzles  0,1,3  and  4  was  measured  from  enlarged 
schlieren  photographs  taken  at  increasing  NPRs  with  no 
secondary  flow.  This  distance  to  the  point  of  flow 
separation  is  labeled  x^.  Nozzle  2  was  not  considered  since 
its  configuration  with  no  secondary  flow  is  identical  to 
nozzle  1.  Figures  7  through  10  are  composites  of  the 
schlieren  photographs  used  to  measure  the  values  of  Xg. 

They  illustrate  the  nozzle  flow  structure  at  increasing 
values  of  NPR.  The  clear  walls  required  for  the  schlieren 
nozzle  assembly  were  made  of  0.25  inch  plexiglass,  not 
optical  glass.  Therefore,  some  inherent  refraction  of  the 
images  was  evident  when  the  enlarged  photographs  were 
examined,  which  degraded  the  accuracy  of  the  measurements. 

Figure  11  is  a  summary  of  all  x^_  values  measured.  The 
nozzle-normalized  length  to  flow  separation  is  simply  x^/L. 
where  L  is  the  length  from  the  nozzle  throat  to  the  nozzle 
exit  along  the  centerline  axis.  The  two  horizontal  lines 
indicate  the  location  of  the  first  secondary  port  Cnot 
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Schlieren  Photographs 


Fig  7  Nozzle  0  Separation 
CNozzle  pressure  ratios  in  order  from  upper  left  hand  to 
lower  right  hand  photograph:  10.3,  11.6,  13.0,  14.2,  15.5, 
17.8,  19.7,  24.6,  29.6,  34.53 
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Fig  8  Nozzle  1  Separation  Schlieren  Photographs 


CNozzle  pressure  ratios  in  order  from  upper  left  hand  to 
lower  right  hand  photograph:  10.3,  12.0,  14.7,  16.1, 
18.6.  20.8.  26.3.  31.8,  38.23 

1  Q 


Fig  9  Nozzle  3  Separation  Schlieren  Photographs 
CNozzle  pressure  ratios  in  order  from  upper  left  hand  to 
lower  right  hand  photograph:  15.4,  17.2,  19.  0,  21.1, 
22.7,  25.6,  28.4,  35.  5,  42.5,  49.  6} 


Fig  lO  Nozzle  4  Separation  Schlieren  Photographs 
C Nozzle  pressure  ratios  in  order  from  upper  left  hand  to 
lower  right  hand  photograph:  19.7.  21.4.  23.  O.  24.6, 
32.5.  40.3.  48.1.  55.  7} 
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CJTVC  NOZZLE  SEPARATION  TESTS 
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Lengths  in  Test  Nozzles 


applicable  for  nozzle  05  .  The  solid  line  labeled  "Thompson 
Separation  Model"  identifies  separation  points  predicted  by 
Thompson's  model  for  separation  in  nozzles  caused  by  adverse 
pressure  gradients  C17:l-85.  Appendix  B  contains  the 
appropriate  equations  and  sample  calculations.  The  distance 
to  points  of  flow  separation  in  nozzles  1.3,  and  4  with 
secondary  ports  tend  to  be  very  close  to  the  secondary 
ports.  Nozzle  0.  which  has  no  secondary  ports,  has  x^/L 
values  slightly  higher  than  those  nozzles  with  secondary 
ports . 

The  proximity  of  x^/L  for  nozzles  with  secondary  ports 
to  the  beginning  of  these  secondary  ports  indicates  that  the 
separation  may  be  shock -i nduced .  If  the  presence  of  the 
secondary  ports  in  the  flow  do  cause  oblique  shocks  to  form, 
as  the  photographs  indicate,  separation  may  occur  slightly 
before  this  shock  C 2: 2175.  This  is  consistent  with  the 
trends  shown  in  Fig  10,  including  the  divergence  of  the 
measured  values  from  the  Thompson  model  C which  only  deals 
with  separation  due  to  adverse  pressure  gradients}.  This 
divergence  is  consistent  with  Cate’s  conclusions  concerning 
the  Thompson  model  C5:475. 

Nozzle  0  appears  to  have  a  combination  of  separation 
mechanisms.  At  NPRs  up  to  approximately  13,  the  separation 
follows  the  Thompson  model  and  therefore  may  be  caused  by 
adverse  pressure  gradient  separation.  At  higher  NPRs, 
normal  shocks  start  to  appear  in  the  nozzle  cavity  and  the 
separation  appears  to  be  caused  by  these  shocks.  Further 


23 


details  of  shock  structures  are  discussed  later. 


Having  flow  separation  induced  by  the  presence  ol  the 
secondary  ports  is  advantageous  since  then  the  point  of 
separation  will  always  be  known.  Also,  it  ensures  that  the 
secondary  flow  will  be  introduced  into  the  separation  region 
of  the  flow  which  is  important  for  vectored  flow  C8:2). 

The  location  of  the  point  of  separation  may  also  be 
influenced  by  the  flow  turning  requirements  of  the  nozzle. 
The  ratio  of  nozzle  exit  area  to  the  area  at  the  secondary 
port  location  is  1.42  for  all  CJTVC  test  nozzles.  The 
oblique  shocks  C which  causes  separation)  observed  may  be 
required  to  turn  the  flow  towards  the  nozzle  exit  as  Fig  12 
illustrates.  Further  discussion  of  nozzle  shock  structures 
and  flow  turning  requirements  follow. 


Shock  Structures  C no  secondary  in  lection) 

The  internal  shock  structures  are  different  for  each 
nozzle  tested.  By  studying  these  structures,  considerable 
insight  is  gained  concerning  the  primary  flow  in  the 
nozzles.  The  following  paragraphs  summarize  the  major  shock 
and  flow  characteristics  of  each  nozzle. 

At  low  NPR ,  oblique  shocks  can  be  seen  in  the  diverging 
portion  of  nozzle  0,  as  shown  In  Fig  7.  These  oblique 
shocks  may  be  caused  by  the  adverse  pressure  gradient 
separation  C15:54).  As  the  NPR  increased,  either  normal 
shocks  or  Mach  reflections  are  clearly  evident.  These  may 
be  caused  by  either  an  intersection  of  the  oblique  shocks 
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Fig  12  Oblique  Shock  Flow  Turning  in  Nozzle: 


for  the  Mach  reflection  case  CB:481D,  or  by  normal  shock 
formation  due  to  increased  over  expansion  Cl 4:  135}.  Both 
mechanisms  will  provide  the  observed  axial  flow  as  Fig  13 
illustrates.  Oblique  shocks  mach  reflections  turn  the  flow 
towards  the  nozzle  exit.  Normal  shock  formation  will 
thicken  the  boundary  layer  and  increase  the  separation 
region  to  form  a  path  that  keeps  the  primary  flow  axial. 

Mach  line  development  will  include  slip  lines  that  may  also 
define  the  primary  flow  path  C 1 4: 557 , 583D .  Since  these 
shock  structures  are  extremely  complicated  to  analyze,  it  is 
impossible  within  the  scope  of  this  study  to  determine  which 
mechanism  causes  the  observed  flow  in  nozzle  O. 

In  nozzle  1  CFig  85,  oblique  shocks  caused  by  the 
secondary  ports  are  evident  in  all  cases.  These  oblique 
shocks  determine  the  flow  turning  direction  as  previously 
illustrated  in  Fig  12.  Increasing  values  of  the  oblique 
shock  angle,  ft,  indicate  increasing  pressure  rises  across 
the  shock.  Therefore,  any  difference  between  oblique  shock 
angles  for  the  upper  and  lower  nozzles  walls  also  means 
there  is  a  difference  in  static  pressures.  These  pressure 
differences  cause  the  separation  region  behind  the  upper 
shock  to  enlarge  and  the  flow  to  attach  to  the  lower  wall  in 
all  cases.  This  biased  flow  is  consistent  with  the  bistable 
nature  of  this  nozzle  configuration,  which  will  be  discussed 
1 ater . 

Nozzle  3  CFig  9}  displays  the  secondary  port-induced 
oblique  shock  structure  similar  to  nozzle  1.  but  has  axial 
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Fig  13  Shock  and  Flow  Structures  for  Nozzle  0 


Primary  Flow 
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flow  stability  similar  to  nozzle  0.  Figure  14  illustrates 
why  nozzle  3  is  not  bistable  Part  Ca}  shows  the  nozzle  as 
it  appears  in  Fig  9.  Subsonic  separation  regions  behind  the 
oblique  shocks  provide  a  path  for  the  supersonic  nozzle 
flow.  If  this  supersonic  flow  tries  to  attach  to  either 
nozzle  wall,  the  secondary  port  located  at  the  nozzle 
maximum  area  will  induce  an  oblique  shock.  Video  tape  flow 
visualization  results  Cdiscussed  in  detail  later}  show  that 
these  shocks  form  with  or  without  secondary  flow.  The  shock 
will  include  a  separated  region  before  and  after  it  C2:217} 
that  will  not  allow  the  flow  to  attach.  Therefore,  the 
presence  of  the  additional  secondary  ports  cause  nozzle  3  to 
have  stable  axial  flow. 

Nozzle  4  CFig  10}  shows  the  same  oblique  shock  structure 
in  the  diverging  section  as  nozzles  1  and  3.  It  has  stable 
axial  flow  like  nozzle  3.  but  for  different  reasons.  Since 
this  nozzle  does  not  have  the  reconverging  section  prior  to 
its  exit,  any  separation  region  caused  by  shocks  will  not  be 
entrained  as  it  is  in  the  CJTVC  design.  At  higher  NPRs , 
small  separation  regions  appear  just  inside  the  nozzle  exit. 
These  may  be  caused  by  the  reflection  of  the  secondary 
port-induced  oblique  shocks  as  illustrated  in  Fig  1?.  At 
lower  NPRs,  the  shock  ana1 e  may  be  of  sufficient  magnitude 
to  cause  the  shock  to  intersect  the  separation  region.  At 
higher  NPRs,  the  shock  angle  is  decreased  and  the  shock 
reflection  is  possible. 
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C a}  Low  NPR  case 


Separati on 
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C  fcO  High  NPR  case 


Fig  15  Shock  and  Flow  Structures  for  Nozzle  4 


Flow  Visualization  Observations 


In  addition  to  the  separation  schlieren  photographs, 
video  tapes  of  schlieren  photography  were  recorded. 
Approximately  two  hours  of  video  tape  was  recorded  in  the 
AFIT  blowdown  wind  turnel  which  allowed  testing  at  reduced 
ambient  pressure  using  the  prototype  nozzle  at  NPRs  between 
4  and  23.  It  had  one  nozzle  half  with  one  port  that 
directed  secondary  flow  perpendicular  to  the  primary  flow 
and  one  half  with  a  port  that  directed  secondary  flow 
perpendicular  to  the  nozzle  wall.  An  additional  40  minutes 
of  video  tape  was  recorded  at  the  VKI  S-l  tunnel  using 
nozzles  0,1,2,  and  3  with  NPRs  between  33  and  41.  This  was 
accomplished  by  focusing  a  video  camera  capable  of  macro 
photography  on  the  schlieren  i mage  formed  on  a  frosted  glass 
plate.  Several  important  observations  are  summarized  below. 

The  bistable  nature  of  nozzles  1  and  2  was  evident  from 
the  video  tapes.  When  primary  flow  was  introduced  into 
either  nozzle,  the  flow  immediately  became  biased  towards 
one  wall.  This  bias  was  more  pronounced  at  lower  NPRs  where 
the  flow  seemed  to  be  more  unstable.  When  secondary  flow 
adequate  for  vectoring  was  introduced,  the  flow  immediately 
attached  to  the  wall  opposite  to  the  operating  secondary 
port.  In  order  to  produce  axial  flow,  both  secondary  flows 
were  required. 

Figure  16  illustrates  shock  structures  found  in  the 
video  tapes  for  vectored  and  axial  flow  cases.  In  the 
vectored  case,  the  secondary  flow  raises  the  pressure  behind 
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its  oblique  shock  and  thereby  increases  the  shock  angle. 

This  increase  in  pressure  on  one  nozzle  half  causes  flow  to 
attach  to  the  other  nozzle  half.  At  the  same  time,  the 
secondary  flow  adds  mass  to  the  separation  region  which  is 
entrained  by  the  shear  layer  caused  by  the  oblique  shock. 

The  combination  of  these  effects  create  a  path  for  the 
nozzle  flow.  For  the  case  of  axial  flow,  a  Mach  reflection 
is  observed  at  the  secondary  flow  location  when  both 
secondary  ports  have  equal  flow.  This  is  due  to  the  oblique 
shock  angles  being  increased  to  where  their  intersection 
forms  a  Mach  reflection  to  satisfy  flow  conditions 
C20:  392, 3933  .  The  combination  of  slip  line  development  and 

increased  regions  of  separated  flow  behind  the  shocks  due  to 
secondary  flow  create  a  path  for  the  nozzle  axial  flow. 
Static  photos  of  this  observation  are  not  available  due  to 
malfunctions  in  the  schlieren  apparatus. 

An  interesting  transient  phenomenon  was  recorded 
accidentally  when  a  failure  in  the  primary  flow  supply  dryer 
allowed  some  water  to  enter  the  nozzle.  As  the  nozzle 
primary  flow  began,  water  droplets  traced  a  vortex  path 
opposite  the  wall  where  the  flow  attached.  Figure  17  Ca3  is 
a  sketch  of  the  path  which  was  observed  on  several 
occasions.  A  possible  mechanism  for  the  vortex  formation  is 
illustrated  in  Fig  17  C  b3  andCc3.  Two-dimensional  vortex 
pairs  have  been  proven  to  form  when  fluid  is  impulsively 
started  through  the  sharp  edges  of  a  channel  opening 
Cl: 1113.  What  occurs  next  is  not  clear;  either  the  flow 
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C b}  Vortex  pair 
generated  by 
initial  flow 


CcD  Vortex  growth  and 


becomes  biased  to  one  wall  causing  the  vortex  opposite  it  to 
increase,  or  some  inequality  between  the  vortex  pair  causes 
the  flow  bias.  A  slight  misalignment  of  the  nozzle  halves 
may  cause  one  vortex  to  begin  before  the  other,  creating 
such  a  bias.  The  entire  start  up  mechanism  of  the  nozzle 
was  too  fast  to  determine  if  the  observed  vortex  trace  was 
the  cause  or  the  effect  of  the  flow  bias.  Future  studies 
using  a  high  speed  camera  may  answer  this  question. 
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IV. 


Nozzle  Performance  Results 


Performance  with  No  Secondary  Flow 

In  order  to  properly  evaluate  nozzle  vectoring 
performance,  it  is  important  to  evaluate  its  performance 
with  only  primary  flow  first.  From  this,  any  inherent 
angular  bias  in  the  flow  can  be  determined  and  taken  into 
account  when  evaluating  vectored  performance.  Also,  effects 
of  vectored  flow  on  nozzle  efficiencies  can  be  evaluated. 

The  VKI  6-component  force  balance  was  used  to  measure 
nozzle  side  and  axial  forces.  From  these  values,  total 
force  and  vector  angles  were  calculated  using  simple 
trigonometric  relations.  Test  NPRs  were  obtained  by  varying 
primary  pressure  between  1.2  and  3.5  bar  absolute  in  a 
reduced  back  pressure  which  varied  between  0.07  and  0.08  bar 
absol  ute. 

Thrust  Force  and  Vector  Angle 

Figures  18  through  22  summarize  the  force  and  vector 

angle  data  for  nozzles  0,1, 3, 4,  and  5.  Nozzle  2  was  not 

tested  since  its  configuration  with  no  secondary  flow  is 

identical  to  nozzle  1.  Each  figure  has  two  sets  of  data 

which  represent  the  results  of  two  separate  test  runs. 

Nozzles  0,1,  and  3  have  approximately  the  same  values  of 
measured  axial  and  total  forces.  Nozzle  4  values  are 
slightly  higher  and  nozzle  5  has  the  highest  values  for 
measured  thrust.  An  examination  of  the  static  pressures 
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VECTOR  ANCLE  (DEC)  MASKED  FORCE  (LBF) 


ozzie 


orces  ana  vector 


C no  secondary  flow} 


VECTOR  ANCLE  (DEC)  KCASURED  FORCE  (LBF) 


C  no  secondary  flowD 
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VECTOR  ANGLE  (DEG)  KASURED  FORCE  (LBF) 


NOZZLE  3  FORCE  MEASUREMENTS 


NO  SECONDARY  FLOW 


Fig  20  Nozzle  3  Forces  and  Vector  Angles 


Cno  secondary  flow} 
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VECTOR  ANGLE  (DEG)  hCASURED  FORCE  (LBF) 


NOZZLE  4  FORCE  MEASUREMENTS 


NO  SECONDARY  FLOW 


10  20  30  40  50  60 


NOZZLE  PRESSURE  RATO  (NPR) 
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inside  each  nozzle  provides  an  explanation  for  the 
variations  in  measured  force  values.  Figure  23  is  a  summary- 
plot  of  the  centerline  static  pressure  values  for  ail 
nozzles.  These  values  were  measured  using  a  scani valve 
connected  to  an  ENDEVCO  pressure  transducer  at  the  same 
conditions  as  the  force  measurements.  Appendix  D  contains 
the  individual  pressure  profiles  for  each  nozzle  used  to 
create  this  summary.  The  average  NPR  of  the  tests  is  42.4  ± 
1 . 2. 

The  pressure  profiles  for  nozzles  0,1,  and  3  all  show  an 
abrupt  pressure  rise  past  the  secondary  port.  This 
indicates  that  a  shock  or  series  of  shocks  are  present.  The 
pressure  rise  for  nozzle  0  is  greater  than  those  for  nozzles 
1  and  3.  The  existence  of  normal  shocks  in  nozzle  0  Cwith 
no  secondary  ports}  and  the  existence  of  oblique  shocks  in 
nozzles  1  and  3  Cwith  secondary  ports}  is  consistent  with 
the  shock  structures  discussed  in  section  III.  This  will 
also  account  for  the  pressure  profiles  since  the  static 
pressure  rise  across  a  normal  shock  is  greater  than  the  rise 
across  an  oblique  shock  for  a  given  initial  Mach  C6:5T8}. 

If  the  normal  and  oblique  shocks  caused  the  only 
pressure  rises  in  the  nozzles,  one  would  expect  nozzles  1 
and  3  to  have  greater  measured  forces  than  nozzle  0  since 
their  total  pressure  loss  is  less.  However,  the  pressure 
profiles  show  additional  pressure  rises  through  the  nozzle 
cavity  for  nozzles  1  and  3.  The  end  result  seems  to  be  that 
the  overall  pressure  losses  in  these  nozzles  are 
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SUMMARY  OF  NOZZLE  CENTERLINE  PRESSURES 


Fig  23  Summary  of  Nozzle  Centerline  Static 
Pressures  C  Aver  age  NPR  =  42.4  ±1.2) 


approximately  the  same. 

The  nozzle  4  pressure  profile  shows  a  rapid  drop  in 
pressure  up  to  the  secondary  port  and  a  subsequent  gradual 
rise  to  the  exit.  This  pressure  drop  indicates  that  the 
waves  previously  observed  in  Fig  10  were  most  likely 
expansion  fans  rather  than  shocks.  This  is  consistent  with 
the  profile  one  nught  expect  from  and  over -expanded , 
converging-diverging  nozzle  C7:160}. 

The  profile  for  nozzle  5  shows  a  smooth  decline  in 
pressure  with  no  pressure  rises  throughout  the  nozzle 
length.  This  indicates  that  no  shocks  were  present  in 
nozzle  5,  and  that  it  was  not  over -expanded .  With  no  total 
pressure  loss  due  to  shocks,  it  is  reasonable  that  nozzle  5 
has  the  largest  measured  force. 

The  vector  angles  of  all  nozzles  are  very  stable  at  NPRs 
greater  then  30.  While  the  angles  are  stable,  they  are  also 
slightly  offset  from  zero  degrees,  which  is  perfectly  axial 
flow.  Several  possible  causes  of  this  slight  offset  for 
each  nozzle  are:  Cl}  misalignment  between  nozzle  halves, 

C2}  manufacturing  differences  between  nozzles  halves,  or  C 3} 
misalignment  between  the  nozzle  assembly  and  the  force 
balance.  Unfortunately,  time  did  not  allow  for  further 
investigation  of  these  offsets. 

Thrust  Efficiencies 

Thrust  efficiency  for  each  nozzle  condition  was 
determined  by  dividing  the  measured  total  thrust  by  the 


ideal  thrust.  Ideal  thrust  calculations  assumed  irentropic 
flow  expansion  where  the  static  pressure  at  the  nozzle  exit 
equals  the  ambient  pressure.  Also,  ideal  exit  velocity  was 
calculated  using  isentropic  relations.  Appendix  B  contains 
complete  details  of  these  equations  and  sample  calculations. 

Figures  24  through  28  show  the  thrust  efficiencies  at 
various  NPRs  for  nozzles  0.1, 3. 4,  and  5.  The  measured  total 
thrust  values  used  to  determine  these  are  from  the  same  data 
as  Figures  18  through  22.  Values  for  nozzles  O  and  3  are 
approximately  the  same,  varying  between  0.66  and  0.78. 
Efficiencies  for  nozzle  1  are  slightly  higher,  varying 
between  0  74  and  O. 87.  There  does  not  appear  to  be  any 
obvious  trend  for  efficiency  versus  NPR  for  these  nozzles. 
Table  1  summarizes  the  maximum  and  minimum  efficiency  values 
for  each  nozzle  along  with  its  average  efficiency  and 
appropriate  standard  deviation.  Nozzle  1  has  the  best 
average  efficiency  of  all  the  nozzles  with  reconverging 
sections . 

Nozzle  4  efficiencies  tend  to  increase  with  increasing 
NPR.  This  is  typical  for  convergi  ng-di  ver  gi  ng  nozzle 
behavior  where  at  lower  NPRs,  the  formation  of  shocks  within 
the  nozzle  may  occur.  As  the  NPR  increases,  the  shock  is 
moved  out  of  the  nozzle  and  the  efficiency  increases  due  to 
the  decrease  in  shock  losses  C6:160D.  Because  of  this  trend 
of  increasing  efficiency,  the  average  efficiency  of  nozzle  4 
is  not  comparable  to  the  more  consistent  values  of  nozzles 
0,1,  and  3.  However,  its  magnitude  indicates  that  there  is 
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m) ST  EFFCENCY  TH?UST  EFFCtNCY 


Fig  24  Nozzle  O  Thrust  Efficiencies 
C  no  secon4a»-y  flow} _ 

NOZZLE  1  THRUST  EFFICIENCIES 
NO  SECONDARY  FLOW 


NOZZLE  PRESSURE  RATIO  (NPR) 


Fig  25  Nozzle  1  Thrust  Efficiencies 
Cno  secondary  flow} 
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ThRUST  EFFCENCY 


NOZZLE  3  THRUST  EFFICIENCIES 
NO  SECONDARY  FLOW 


Fig  26  Nozzle  3  Thrust  Efficiencies 
Cno  secondary  flow} 


NOZZLE  4  THRUST  EFFICIENCIES 
NO  SECONDARY  FLOW 


NOZZLE  PRESSURE  RATO  (NPR) 

Fig  27  Nozzle  4  Thrust  Efficiencies 

Cno  secondary  flow} 
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THRUST  EFFCENCY 


Fig  28  Nozzle  5  Thrust  Efficiencies 
Cno  secondary  flow) 

Table  1,  Nozzle  Thrust  Efficiencies  with  No  Secondary  Flow 


Nozzle 


Mi  ni  mum 

Maxi  mum 

Aver  age 

Ef  f  i  ci  ency 

Ef f i ci ency 

Ef  f i ci ency 

0.  664 

0.  763 

O.  72510-  025 

0.  737 

0.  870 

O.  76510.  034 

0.  664 

0.  778 

O.  73310.  034 

0.  625 

0.  81  3 

0.  70610.  058 

0.  852 

0.  948 

0.  91310.  026 

little  efficiency  penalty  for  adding  the  reconverging 
section  to  the  convergi ng-di vergi ng  design  of  nozzle  4. 
Nozzle  5  efficiencies  are  considerably  greater  than  all 
other  nozzles.  Pressure  data  for  nozzle  5  indicates  that  it 
does  not  suffer  from  shock  losses. 

Performance  with  Secondary  Flow 

Unless  otherwise  noted,  all  tests  with  secondary  flow 
were  performed  with  an  applied  primary  pressure  of  2.923  bar 
absolute  at  the  nozzle  entrance.  This  was  the  largest 
pressure  that  could  be  supplied  consistently  from  the 
existing  test  facilities  to  the  nozzle.  The  operating  NPR 
varied  slightly  according  to  variations  in  the  back 
pressure.  Tests  were  performed  by  keeping  the  primary 
pressure  constant  and  varying  the  mass  flow  ratio  CMFRD  of 
secondary  to  primary  mass  flow  rates.  Appendix  C  provides 
the  pressure  and  mass  flow  values  for  all  test  conditions. 
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1ESS  SDE  FORCE  (Fs/PpAt) 
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DMEN30NLESS  SDE  FORCE  (Fs/PtAl) 
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C*€WSCMESS  SCE  FORCE  (Fs/PtAl) 


Fig  31  Nozzle  3  Vector  Angle  and  Side  Force  Variations 


ESS  SCE  FORCE  (Fs/PtAt)  VECTOR  ANCLE  (DEC) 


Fig  32  Nozzle  4  Vector  Angle  and  Side  Force  Variations 


ports,  “+1  SI"  is  the  upper  wall  secondary  port  and  "-1  SI" 
is  the  lower  wall  port. 

All  vector  angle  plots  show  positive  angle  bias,  that 
is,  greater  values  of  positive  vector  angles  than  negative 
vector  angles.  Some  of  this  may  be  due  to  misalignment 
between  the  nozzle  assembly  and  force  balance,  but  the 
consistency  of  the  bias  direction  indicates  that  it  is  more 
likely  due  to  some  slight  variation  in  the  nozzle  test 
assembly  itself. 

Figure  29  shows  that  the  nozzle  1  produces  vectored 
thrust  at  all  MFRs  tested.  The  vector  angles  produced  with 
secondary  flow  from  one  wall  decrease  slightly  with 
increasing  MFR.  The  axial  flow  produced  using  secondary 
flow  from  both  walls  has  slightly  less  positive  offset  from 
zero  degrees  than  the  nozzle  had  with  no  secondary  flow  CFig 
193. 

Figures  30  and  31  show  that  nozzle  2  and  3  have  definite 
transition  points  between  unvectored  and  vectored  flow. 
However,  these  transition  points  do  not  occur  at  the  same 
MFR  for  both  vector  directions,  again  indicating  that  some 
internal  variations  in  the  nozzles  may  account  for  this. 

When  vectored  flow  is  achieved,  the  vector  angles  produced 
remain  essentially  constant  with  increasing  MFR.  As  with 
nozzle  1,  the  axial  flov,-  produced  using  secondary  flow  from 
both  walls  for  nozzles  2  and  3  has  less  positive  offset  from 
zero  degrees  than  the  nozzles  had  with  no  secondary  flow 
CFigs  19  and  203. 
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Figure  32  shows  that  nozzle  4  produced  relatively  small 
vector  angles,  even  at  high  MFR.  This  illustrates  the 
ef f ecti veness  of  the  reconverging  section  in  producing 
vectored  flow. 

For  all  figures,  the  dimensionless  side  force  values 
seem  to  level  off  as  the  secondary  to  primary  pressure  ratio 
increases.  Nozzles  2  and  3  in  particular  show  sharp  side 
force  value  rise  to  a  plateau,  similar  to  trends  found  by 
Lambert  and  Franke  Cll:43.  Figure  33  shows  the  relation 
between  MFR  and  the  secondary  to  primary  pressure  ratios 
used  in  Figs  29  through  32. 

Table  2  summarizes  the  average  vector  angle  resulting 
from  secondary  flow  from  one  wall  measured  from  the  angle 
resulting  from  secondary  flows  from  both  walls.  In  most 
cases,  the  vector  angle  is  biased  to  the  positive  direction. 

Figure  34  is  a  summary  of  all  data  contained  in  Figs  29 
through  32.  It  presents  the  data  as  vector  ranges,  which 
are  simply  the  arithmetic  difference  between  the  two  vector 
angles.  Nozzle  1  has  the  best  vector  range  values  at  MFRs 
below  0.147  while  nozzle  2  has  the  best  vector  range  values 
above  that  MFR.  While  nozzles  2,3,  and  4  all  have 
increasing  vector  ranges  with  increasing  MFR,  the  range  of 
nozzle  1  vector  angles  decreases  slightly  with  increasing 
MFR.  The  trend  of  nozzle  1  is  consistent  with  the  trends  of 
Talda’s  data  with  a  similar  nozzle  C16:58-6CD. 

The  MFR  values  which  produced  the  best  vector  ranges  at 
the  flow  conditions  of  a  primary  pressure  of  2.923  bar 
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SECONDARY: PRM ARY  MASS  FLOW  RATIO 


MASS  FLOW  RATIO  VERSUS  PRESSURE  RATIO 


(SECOND ARY:PRIM ARY  FLOWS) 


SECONDARY:PRWARY  PRESSURE  RATO 
Fig  33  Mass  Flow  Ratio  versus  Pressure  Ratio 
C  Secondary  Flow  to  Primary  Flow} 

Table  2.  Vector  Angle  Measured  from  Angle  Produced  by 

Both  Secondary  Flows 


MFR 

NOZZLE  1 

NOZZLE  2 

NOZZLE  3 

NOZZLE  4 

0.  067 

+  14.  9/-18.  2 

+  15.  5/+2.  5 

+1 . 7/-1 . 2 

+1  .  1/+0.  1 

0.  115 

+  14. 7/--17.  5 
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+ 

v-» 

\ 

1 

O 

(D 

— 

0.  147 
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— 

0.  239 
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VECTOR  ANCLE  (DEG) 


VECTOR  ANGLE  RANGES  FOR  ALL  NOZZLES 
(TESTED  AT  Pp  =  2.92  BAR  ABS) 
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Fig  35  Vector  Range  Variation  with  NPR 
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absolute  were  selected  for  further  testing.  For  this 
testing,  the  MFR  value  for  the  appropriate  nozzle  remained 
constant  while  the  NPR  varied.  Table  3  lists  the  MFR  values 
and  their  resulting  best  vector  ranges  for  all  nozzles. 


Table  3.  Best  Vector  Angle  Range  Values 
for  Nozzles  Tested 


Nozzl e 

MFR 

Best  Vector  Range 

1 

0.  087 

33.  1 

2 

0.  193 

29.  5 

3 

0.  291 

18.  3 

4 

0.  239 

7.  4 

Figure  35  summarizes  the  vector  angle  range  variance 
with  NPR.  With  varying  NPR,  all  nozzles  keep  their  relative 
position,  nozzle  1  had  the  best  vector  range,  followed  by 
nozzle  2  and  then  nozzle  3.  All  nozzles  show  little  range 
variation  with  varying  NPR,  except  nozzle  3  which  has 
decreased  range  with  increasing  NPR. 

Thrust  Efficiencies 

Thrust  efficiencies  for  flow  with  secondary  injection 
were  determined  by  the  same  methods  used  for  flow  with  no 
secondary  injection,  except  the  ideal  thrust  was  calculated 
by  multiplying  the  isentropic  exit  velocity  by  the  total 
mass  flow  rate  Cprimary  plus  secondary).  Figures  36  through 
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39  show  the  thrust  efficiencies  at  various  MFRs  for  nozzles 


1,2,3,  and  4.  In  some  cases,  the  nozzles  show  a  slight 
decrease  in  efficiency  with  increasing  MFR.  This  decrease 
may  be  caused  Dy  increased  nozzle  losses  as  Increasing  MFR 
is  accompanied  by  an  increased  ratio  between  secondary  and 
primary  pressures.  This  increased  pressure  ratio  can  lead 
to  increased  shock  strengths  and  their  resulting  losses  in 
the  nozzle.  In  most  cases,  the  thrust  efficiencies  varied 
less  than  through  the  MFR  range  tested. 

Table  4  shows  the  average  thrust  efficiencies  for  all 
nozzles  at  the  three  vectoring  conditions.  Nozzle  1  had  the 
best  efficiency  with  the  no  secondary  flow  case,  and  it  has 
the  best  efficiencies  for  individual  secondary  flow 
injection  cases.  Nozzle  3.  however,  has  the  best  efficiency 
when  both  secondary  flow  are  used.  The  average  efficiencies 
determined  for  these  CJTVC  nozzles  are  in  excellent 
agreement  with  the  efficiencies  determined  by  Talda  for 
vectored  flow  also  using  a  circular  secondary  injection  port 
C  16:  973  . 

Side  Momentum  Gains 

As  previously  mentioned,  the  side  momentum  gain  CG^3  is 
defined  as  the  ratio  of  the  measured  side  force  to  the 
momentum  of  the  secondary  jet  flow.  Obviously,  large  values 

of  G  are  desirable  since  it  reflects  that  a  large  side 

M 

force  is  being  produced  by  a  small  quantity  secondary  jet 
momentum.  Since  the  side  force  is  a  result  of  primary  and 
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VECTORED  T>«UST 


Fig  36  Nozzle  1  Thrust  Efficiencies 


Cwith  secondary  flow? 


Fig  37  Nozzle  2  Thrust  Efficiencies 


Cwith  secondary  flow} 
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VECTORED  TVRUST  0TEENCY  VECTORED  TVRLKT 


Fig  38  Nozzle  3  Thrust  Efficiencies 


Cwith  secondary  flow} 


Fig  39  Nozzle  4  Thrust  Efficiencies 


Cwith  secondary  flow} 
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Table  4.  Nozzle  Thrust  Efficiencies  with  Secondary  Flow 


NOZZLE 

+  1  SI 

-1  SI 

BOTH  SI 

1 

0.  780  ±  0.  013 

C.  778  ±  0.  0006 

0.  774  ±  0.  0009 

2 

O.  767  ±  O.  012 

O.  769  ±  0.  Ol  9 

0.  770  ±  O.  018 

3 

O.  776  ±  O.  009 

O.  771  ±  0.  OlO 

0.  793  ±  O.  016 

4 

O.  742  ±  O.  013 

0.  744  ±  0.  017 

0.  753  ±  O.  012 
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secondary  flow  interactions,  it  is  also  of  interest  to 

consider  the  ratio  of  the  secondary  flow  momentum  to  the 

primary  flow  momentum  CR  ).  Appendix  B  includes  the 

M 

equations  and  sample  calculations  for  determining  values  of 
GM  and  RM 

Figures  40  through  43  show  the  values  of  Gu  and  R  at 

M  M 

varying  MFR  for  nozzles  1,2,3,  and  4.  R^  increases  linearly 
with  increasing  MFR,  indicating  that  the  secondary  velocity 
is  approximately  constant. 

If  the  vectored  flow  were  dependent  only  on  the  momentum 

of  the  secondary  flow,  one  would  expect  that  G..  would 

M 

increase  with  increasing  R  However,  the  opposite  is  true 
for  nozzles  1  and  2.  indicating  that  something  more  than  the 
"brute  force"  mechanism  applies.  The  trends  in  Figs  40  and 
41  agree  with  Talda’s  data  that  shows  side  force  increasing 
to  a  maximum  and  then  decreasing  with  further  increases  in 
MFR  Cl  6: 54-56J) . 

Nozzles  3  and  4  show  very  small  values  of  G  which 
remain  essentially  constant  with  increasing  MFR.  This  is 
consistent  with  the  poor  vectoring  performance  of  these 
nozzles.  As  discussed  in  section  III,  the  existence  of  an 
additional  secondary  flow  port  on  the  curved  wall  of  nozzle 
3  may  prevent  the  supersonic  vectored  flow  from  attaching, 
thereby  leading  to  low  G  The  use  of  suction  at  this 
secondary  port  may  allow  for  flow  attachment,  but  that  is 
beyond  the  scope  of  this  study. 


63 


SECtfCAf^MAtfY  MOMENTUM  RATO 


SECONDARY PRMARY  MOMENTUM  RATE)  SDE  MCMJTFIM  GAN 
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SECO€ARYfRMARY  MOMENTUM  RATO 
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SECOfCARYPHUARY  VXMMTUM  RATD 


NOZZLE  4  FLOW  MOMENTUM  RATD 
AVERAGE  NPR  =  37.0 


Fig  43  Nozzle  4  Momentum  Gains  and  Flow 
Momentum  Ratios 
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Static  Pressure  Profiles 


Figures  44  through  46  show  the  centerline  static 
pressures  for  nozzles  1,2,3,  and  4  measured  at  the  best 
vector  range  condition  for  each  nozzle  as  listed  in  Table  3. 
These  figures  can  be  compared  with  Fig  23  which  provides  the 
centerline  static  pressure  profile  for  nozzles  with  no 
secondary  flow. 

Figures  44  and  45  are  the  pressure  profiles  for 
secondary  flow  from  only  one  nozzle  wall.  As  expected, 
these  two  figures  are  almost  identical.  The  vectored  flow 
profiles  for  nozzle  1  differ  from  its  profile  with  no 
secondary  flow,  as  illustrated  in  Fig  23.  Both  profiles 
show  shar  p  pressure  rise  close  to  the  secondary  flow  port, 
but  the  pressure  rise  for  the  vectored  case  occurs  closer  to 
the  throat.  The  pressure  values  after  the  sharp  rise  drop 
slightly  and  then  remain  essentially  constant  until  an 
expansion  Is  required  at  the  exit.  The  difference  in  the 
pressure  rise  location  between  the  cases  with  and  without 
secondary  flow  is  consistent  with  the  mechanism  previously 
illustrated  in  Fig  16.  The  introduction  of  secondary  flow 
raises  the  static  pressure  in  the  area  behind  the  oblique 
shocks  in  the  nozzle  diverging  section.  This  change  in  flow 
conditions  causes  the  oblique  shock  angle  to  increase,  which 
in  turn  causes  the  centerline  pressure  rise  due  to  the  shock 
to  move  closer  to  the  throat,  as  observed.  The  nozzle  2 
static  pressure  profile  is  similar  to  the  nozzle  1  profile, 
except  it  has  slightly  higher  values.  The  higher  values 
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reflect  the  higher  secondary  flow  pressures  required  for  the 
maximum  vector  range  for  nozzle  2. 

For  the  case  with  no  secondary  flow,  the  profile  for 
nozzle  3  was  essentially  the  same  as  that  for  nozzle  1. 

With  secondary  flow,  the  profiles  are  different  in  that 
nozzle  3  has  a  more  gradual  pressure  rise  in  the  diverging 
section  than  nozzles  1  and  2.  The  nozzle  3  pressure  profile 
also  has  slightly  higher  values  than  that  for  nozzle  2, 
which  reflects  its  higher  secondary  flow  pressure 
requirements  for  maximum  vector  range. 

The  pressure  profiles  with  secondary  flow  for  nozzle  4 
show  abrupt  pressure  rises  similar  to  the  rises  in  nozzles  1 
and  2.  These  rises  occur  immediately  downstream  of  the 
secondary  flow  ports.  Since  a  similar  pressure  rise  did  not 
occur  in  the  profile  with  no  secondary  flow,  it  appears  that 
the  introduction  of  secondary  flow  into  nozzle  4  produces  a 
shock . 

Figure  46  shows  the  pressure  profiles  for  when  secondary 
flow  is  applied  from  both  nozzle  walls.  The  profiles  with 
both  secondary  flows  are  similar  to  the  profiles  with  single 
side  secondary  flows  for  all  nozzles  at  x/L  values  between 
0.3  and  1.0.  The  profiles  for  both  secondary  flows  have 
higher  pressure  values,  consistent  with  the  resulting  higher 
mass  flow.  Some  changes  are  evident  in  the  area  around  the 
secondary  injection  ports.  When  secondary  flows  from  both 
walls  are  introduced,  all  nozzles  have  abrupt  pressure  rises 
which  occur  closer  to  the  nozzle  throat  than  they  did  for 
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the  cases  of  secondary  flow  being  introduced  from  only  one 
wall.  This  is  consistent  with  the  formation  of  normal  shock 
or  Mach  reflection  in  the  diverging  portion  of  the  nozzle 
previously  illustrated  in  Fig  16. 

In  summary,  the  centerline  static  pressure  profiles  for 
cases  of  secondary  flow  from  either  one  or  both  nozzle  walls 
confirm  the  previously  discussed  flow  visualization 
obser  vati ons . 

Nozzle  Switching  Response 

The  force  and  vectoring  angle  performance  of  CJTVC 
nozzles  have  been  discussed.  If  a  thrust  vectoring  nozzle 
is  to  be  used  as  some  form  of  vehicle  guidance,  its 
switching  response  time  must  also  be  known.  The  switching 
response  time  is  the  time  required  f or  a  change  in  secondary 
flow  to  create  a  change  in  the  nozzle  exit  flow.  These 
times  were  measured  using  a  TEKTRONIX  model  3430  digital 
oscilloscope,  which  simultaneously  recorded  voltage  signals 
from  both  side  force  channels  of  the  VKI  6-component  force 
balance.  The  oscilloscope  started  recording  data  when  an 
outside  trigger  signal  was  applied.  This  trigger  signal  was 
provided  through  the  secondary  flow  solenoid  valve  controls. 
After  the  signal  is  applied  and  the  data  is  recorded,  the 
oscilloscope  digitizes  and  displays  the  signal  voltage 
versus  time  plot  for  analysis.  Cursor  movement  on  the 
oscilloscope  provides  the  capability  for  accurate  time 
measurements . 
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Although  the  oscilloscope  t rigger  signal  and  the 

solenoid  open  or  close  command  signal  were  simultaneously 

applied,  some  delay  time  existed  until  the  secondary  flow 

started  or  stopped  flowing  in  the  nozzle  assembly.  The 

delay  time  consists  of  the  time  required  for  the  solenoid 

valve  to  open  or  close  once  it  received  the  command  and  the 

time  required  for  the  appropriate  flow  change  to  occur  in 

the  secondary  flow  supply  lines  running  to  the  secondary 

ports  in  the  nozzle  assembly.  This  delay  time  was  measured 

by  placing  an  ENDEVCO  pressure  transducer  in  the  secondary 

supply  line  just  prior  to  its  entry  into  the  nozzle 

assembly.  The  output  signal  of  the  transducer  was  recorded 

by  the  oscilloscope  when  the  appropriate  solenoid  command 

was  given.  Figure  47  shows  the  oscilloscope  traces  recorded 

by  a  pen  plotter  for  solenoid  open  and  close  commands.  Two 

times  were  measured  for  each  case.  The  first  time  Ct.3  was 

A 

the  delay  time  for  the  solenoid  valve  to  open  or  close, 
which  was  measured  at  the  point  where  the  transducer  voltage 
slope  started  to  change.  The  second  time  Ct_p  was  the  delay 

D 

time  for  the  secondary  pressure  to  reach  steady  state.  It 

was  measured  at  the  port  where  the  transducer  voltage  slope 

leveled  off.  For  the  valve  open  command,  t  =  25.9  ms  and 

tg  =  37.8  ms,  while  the  valve  close  command  had  t^  =  37.8  ms 

and  t  =  79. 8  ms. 
ti 

Switching  response  times  were  measured  from  the  strain 
gage  voltage  signals  of  the  two  side  force  channels  of  the 
force  balance.  Appendix  A  contains  information  concerning 
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Ca?  Valve  open  command 
del  ays 

t .  =  25. 9  ms 

A 

1 0  =  65.  9  ms 
B 


-  delay  time  for  solenoid  valve  to  open  or  close 
A 

t  =  delay  time  for  secondary  supply  line  pressure  to 
B 

increase  or  decrease  to  steady  state 


Cb)  Valve  close 

command  delays 

t  =37.8  msec 
A 

t.  =  79.  8  msec 
B 


Fig  47  Secondary  Flow  Response  Time  Delays 
(y-axis  divisions  =  50  mV;  x-axis  divisions  =  20  ms!? 
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the  location  of  these  strain  gages.  Time  restrictions  only 
allowed  one  nozzle  configuration  to  be  tested.  Nozzle  1  was 
selected  since  it  had  the  best  vectored  thrust  efficiencies 
and  momentum  gains  of  the  CJTVC  nozzles.  Response  times  for 
six  different  switching  cases  were  measured.  The  secondary 
flow  conditions  before  and  after  the  appropriate  solenoid 
command  for  each  case  are  summarized  in  Table  5.  Response 
times  for  each  case  were  measured  using  MFR  values  of  0.087, 
0.115.  and  0.  147. 

Figures  48  and  49  are  typical  voltage  versus  time  plots 
for  the  six  test  cases.  The  traces  show  the  voltage  for 
each  channel  just  before  the  trigger  signal  is  received  and 
the  resulting  change  in  side  force  voltages  from  the  change 
in  secondary  flow.  When  the  sioe  force  voltages  change 
magnitude,  they  start  to  oscillate  at  a  frequency  of 
approximately  19  Hz.  To  investigate  this  oscillation,  an 
impulse  was  applied  to  the  6-component  balance  with  its 
voltage  output  being  input  to  a  Hewlett-Packard  Model  3582A 
spectrum  analyzer.  The  resulting  power  spectral  density 
plots  showed  a  spike  at  19.2  Hz  after  the  impulse  was 
applied  which  was  not  present  before  the  impulse  was 
applied.  This  was  repeated  several  times  with  the  same 
result,  thereby  confirming  that  the  side  force  oscillations 
in  Figs  48  and  49  are  due  to  the  natural  frequency  of  the 
force  balance. 

Response  times  were  measured  from  the  trigger  signal 
indication  point  to  the  point  where  the  side  force  channel 
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Table  5.  Nozzle  Flow  Vectoring  Response  Time  Test 

Condi ti ons 


Test 

Case 


Flow  Condition 

Before  Solenoid 
Command 

Solenoid  Valve 
Command 

Both  secondary 
flows  applied 

Close  valve  to 

1 ower  wal 1 
secondary  flow 

Beth  secondary 
flows  applied 

Close  valve  to 
upper  wal 1 
secondary  flow 

Lower  wal 1 
secondary  flow 
appl i ed 

Open  valve  to 
upper  wal 1 
secondary  flow 

Upper  wall 
secondary  flow 
appl l ed 

Open  valve  to 
lower  wall 
secondary  flow 

Lower  wall 
secondary  flow 
appl i ed 

Si mul taneousl y 
close  valve  to 
lower  secondary 
flow  and  open 
valve  to  upper 
secondary  flow 

Upper  wall 
secondary  flow 
appl  i  ed 

Si mul taneousl y 
close  valve  to 
upper  secondary 
flow  and  open 
valve  to  lower 
secondary  flow 

Nozzle  Exit  Flow 
Di  recti  on 
1  *  Bef ore  Command 
-  After  Command 
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Fig  48  Nozzle  1  Side  force  Oscilloscope  /races 


Test  Cases  1  -3 

CMFR  =  0.087;  y-axis  d  visions  =  50  mV;  x-axis 
divisions  =  50  ms,  side  force  channel  1  is  top 
lr#c«|  side  force  channel  2  is  bottom  trace} 
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Test  Case  4 
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Test  Case  5 


Test  Case  6 


Fig  49  Nozzle  1  Side  Force  Oscilloscope  Traces 

Test  Cases  4-6 

CMFR  =  0.087,  y-axis  divisions  =  50  mV,  x-axis 
divisions  =  50  msec;  side  force  channel  1  is  top 
trace;  side  force  channel  2  is  bottom  trace) 


traces  begin  a  repeatable  oscillation  pattern.  The 
appropriate  solenoid  valve  command  and  pressure  change  delay 
times  were  subtracted  from  these  times.  The  resulting  times 
represent  Cl)  the  response  time  from  when  the  solenoid  valve 
opens  or  closes  until  the  side  forces  change  and  C  2)  the 
response  time  from  when  the  secondary  flow  reaches  its 
supply  line  pressure  until  the  side  forces  change.  Table  6 
presents  these  response  time  values  for  the  test  cases  and 
MFRs  previously  mentioned.  A  MFR  of  0.  115  produces  the 
shortest  response  times  for  test  cases  1  and  2.  For  test 
cases  3  and  4,  the  response  times  C t^D  produced  by  MFRs  of 
0.  067  and  0.115  are  approximately  equal,  and  are  much  lower 
than  the  times  for  an  MFR  of  O.  147.  Response  times  for  test 
cases  5  and  6  are  shortest  for  and  MFR  of  0.147.  No 
consistent  trends  between  MFR  values  and  response  times  are 
obvious.  The  response  frequencies  Cf^)  for  nozzle  1  at  its 
best  vector  range  MFR  of  0.  067  vary  between  7.  13  and  23.  3  Hz 
when  measured  from  solenoid  valve  opening  or  closing  and 
vary  between  10.2  to  333.3  Hz  when  measured  from  the  point 
when  supply  line  pressure  is  achieved  at  the  secondary  port. 

Summary 

Table  7  summarizes  all  of  the  performance  parameters 
presented  in  this  section  for  test  nozzles  at  the  MFR  that 
produced  the  largest  vector  angles.  As  previously 
discussed,  nozzle  1  is  superior  in  all  aspects  tested  except 
for  thrust  efficiency  with  secondary  flow  from  both  walls. 
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Table  6 


Nozzle  1  Flow  Vectoring  Response  Times 


and  Frequencies 


MFR 

=0.  087 

MFR 

=0.115 

MFR 

=0  1  47 

Test 

Del  ay 

tR 

fR 

tR 

fR 

tR 

fR 

Case 

1 1  me 

C  ms} 

C  Hz} 

C  ms} 

C  Hz} 

C  ms} 

C  Hz} 

1 

tA 

72.  8 

13  7 

69.  5 

PH 

92.  0 

10.9 

S 

30.  8 

32.  5 

27.  5 

36.  4 

70.  O 

14.3 

2 

tA 

112.7 

8.  9 

77.  8 

12.  9 

9.  9 

S 

80.  7 

12.  4 

35.  8 

27.  9 

u 

16.  9 

3 

lA 

43.  O 

23.  3 

45.  5 

22.  0 

87.  5 

11.4 

lB 

3.  0 

333.  3 

5.  5 

181 . 8 

27.  5 

26.  7 

4 

tA 

45.  3 

22.  1 

45.  0 

22.  2 

85.  O 

11.8 

5.  3 

188.  7 

5.  O 

200.  0 

O 

in 

22.  2 

5 

74.  5 

13.  4 

79.  3 

12.  6 

49.  0 

20.  4 

32.  5 

30.  8 

37.  3 

26.  8 

7.  0 

142.  9 

6 

tA 

140.  2 

7.  1 

82.  8 

12.  1 

78.  5 

12.  7 

98.  2 

10.  2 

40.  8 

24.  5 

36.  5 

27.  4 

Table  7  Summary  of  Nozzle  Best  Vector  Range  Performance 
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Nozzl e 

MFR  for 
Best 
Vector 
Range 

Best 

Vector 

Best  Momentum 
Gai  n 

Ef f  i  ci  enci es 
<?  Best  MFR 

Angl  e 
Range 

+1 

-1 

+1 

-1 

2 

1 

O.  087 

33.  1 

4.  50 

3.  35 

0.  794 

0.  784 

O.  786 

2 

O.  193 

29.  5 

1 . 72 

1 . 38 

0.  767 

O.  777 

O.  765 

3 

0.  291 

18.3 

O.  91 

0.  58 

O.  783 

O.  784 

O.  797 

4 

O.  239 

7.  4 

0  5 

0.  1 

0.  737 

O.  740 

0.  753 

+1  =  upper  wall  secondary  flow  applied 
-1  =  lower  wall  secondary  flow  applied 
2  =  both  walls  secondary  flow  applied 
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However,  its  thrust  efficiency  for  this  case  is  only  1.4  j 

percent  less  than  the  best  efficiency  for  the  same  case.  ; 

i 

-I 

Again,  since  the  tested  parameters  for  nozzle  1  showed 
it  to  have  the  best  overall  performance,  its  external  flow 
patterns  were  photographed  using  schlieren  techniques. 

Figure  50  is  a  composite  of  these  photographs  for  MFR  values 
of  O.  087,  O.  115,  0.  147,  and  O.  239.  For  each  MFP,  there  are 

three  photographs,  one  for  each  flow  condition  Csecondary 
flow  from  upper  wall,  secondary  flow  from  both  walls, 
secondary  flow  from  lower  wall  3.  The  horizontal  line  in 
each  photograph  is  a  reference  line  drawn  perpendicular  to 
the  nozzle  exit,  thereby  indicating  the  zero  degree  flow 
angl  e. 
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V. 


Nozzle  Anal  vti cal  Model  Results 

* 

The  previous  sections  of  this  report  have  described  the 
results  of  operating  mechanism  and  performance  tests  This 
section  deals  with  attempts  to  analytically  model  nozzle 
performance  Several  existing  models  are  presented  for 
predicting  the  side  force  to  axial  force  ratio.  The  author 
presents  a  successful  model  for  predicting  the  nozzle  total 
thr  ust 


Viscous  Effects  Analysis 

Before  any  simple  analytical  modeling  is  considered,  it 
is  important  to  determine  if  viscous  effects  in  the  nozzle 
should  be  considered.  Boundary  layer  thickness  will  reach 
its  maximum  value  before  the  flow  separates  from  the  nozzle 
wall  in  its  diverging  section.  Since  previous  discussions 
in  this  report  indicate  that  flow  separation  occurs  near  the 
secondary  flow  ports,  the  value  of  boundary  layer  thickness 
will  be  calculated  at  this  point. 

Both  length  and  hydraulic  di ameter -based  Reynolds  number 
calculations  predict  turbulent  flow  at  the  secondary  points. 
The  walls  in  the  diverging  section  prior  to  the  secondary 
ports  are  flat.  Therefore,  the  following  equation  for  the 
boundary  layer  from  turbulent  flow  over  a  flat  plate  is 


used: 


«5  C.  x} 


0.  37  x 


Cl} 
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where  6^  is  the  boundary  layer  thickness,  x  is  the  length 
along  the  flat  plate,  u^  is  the  flow  velocity,  and  u  is  the 
flow  kinematic  viscosity  C13:633_>.  Using  the  same  model, 
the  displacement  thickness,  6^,  and  the  momentum,  6^,  can  be 
cal  cul  ated: 

<5  =  0.  125  <5  C2D 

<52  =  0.097  6t  C  3D 

where  the  displacement  thickness  is  a  measure  of  the 
displacement  of  the  flow  resulting  from  the  presence  of  the 
flat  plate  and  its  boundary  layer.  The  momentum  thickness 
is  the  measure  of  the  momentum  flux  decrement  caused  by  the 
boundary  layer  ClO:50D. 

Assuming  isentropic  flow  up  to  the  secondary  port,  the 
resulting  values  were  6^  =  0.00548  in,  6^  =  0.00069  in,  and 

6^  =  0.00053  in.  Considering  the  nozzle  dimensions  at  the 
secondary  port,  these  values  show  the  calculated  boundary 
layer  value  is  2.1  percent  of  the  nozzle  height  and  2.2 
percent  of  the  nozzle  width.  Therefore,  it  occupies  8.5 
percent  of  the  nozzle  area.  Since  these  values  are 

relatively  small,  further  consideration  of  the  boundary 
layer  thickness  was  not  included  in  the  analytical  models. 
Detailed  boundary  layer  sample  calculations  are  included  in 
Appendix  E. 

Side  Force  to  Axial  Force  Ratio  Models 

Several  prediction  models  for  side  force  to  axial  force 
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ratios  for  TVC  nozzles  have  been  developed  from  experimental 
results.  Zukoski  and  Spaid  proposed  a  model  which  had  the 
side  to  axial  force  ratio  varying  directly  with  the  ratio  of 
secondary  mass  flow  to  primary  mass  flow  C  20:  1696}.  The 
results  of  section  IV  of  this  report  indicate  that  this 
model  is  not  applicable. 

Broadwell  developed  a  model  using  blast  wave  theory  that 
accounts  for  secondary  flow  momentum  deficit  and  volume 
addition  effects  His  side  to  axial  force  ratio  equation 

l  s : 
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where  cK  is  an  energy  term  defined  by  Broadwell,  M  ,  V  , 

00  CD 

and  1.  are  the  primary  flow  Mach  number,  velocity,  and 
too 

total  temperature  at  the  secondary  port,  y  is  the  specific 

heat  ratio  lor  the  primary  flow,  V  is  the  effective 

ex 

secondary  exhaust  velocity  to  vacuum,  is  the  effective 

secondary  exhaust  velocity,  T  is  the  secondary  flow  total 
temperature,  and  rrv  are  the  primary  and  secondary  mass 

flow  rates,  a  is  the  angle  between  the  secondary  and  primary 
axes,  and  (i  is  the  angle  between  a  line  normal  to  the  nozzle 
wall  at  the  secondary  port  and  the  primary  axis  C3:1071Z). 
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By  considering  the  rates  of  change  of  primary  flow 
momentum  and  secondary  flow  momentum  along  with  experimental 
results.  Green  and  McCullough  developed  a  model  for  fluid 
injection  thrust  vector  control.  Their  side  to  axial  force 
ratio  equation  is. 


P  J 


0.  5 


sp 


0.  5 


l  a  pi  ] 


O.  £5 


C5} 


where  w^  and  w^  are  the  secondary  and  primary  weight  flow 
rates.  I ^  is  specific  impulse,  A  p^  is  the  secondary  flow 
diflerence  between  supply  and  exit,  and  C  is  an  empirical 
constant.  (.9:  577}.  For  constant  mass  flow  and  thrust, 
specific  impulse  is  defined  as: 


I 

sp 


C  6} 


where  F^  is  total  thrust  and  w^  is  total  weight  flow  rate 

Cl  5:  £2} . 

Figures  51  through  54  show  the  relation  between  the 
differences  between  the  measured  side  to  axial  force  ratios 
and  the  ratios  as  modeled  by  the  Broadwell  and 
Green/McCul 1 ough  equations  for  nozzles  1,2,3,  and  4. 
Appendix  B  contains  sample  calculations  for  these  values. 
For  brevity,  "Green"  appears  on  the  graphs  instead  of 
"Green/McCullough".  Clearly,  the  Broadwell  model  does  not 
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SOE  TO  AXIAL  FORCE  RATO  SOE  TO  AXIAL  FORCE  RATOS 


NOZZLE  1  PREDICTION  MODELS 


SOE  TO  AXIAL  FORCE  RATIOS 


Fig  51  Nozzle  1  Side  to  Axaal  Force  Ratio 
Prediction  Results 


NOZZLE  2  PREDICTION  MODELS 
SIDE  TO  AXIAL  FORCE  RATIOS 


Fig  52  Nozzle  2  Side  to  Axial  Force  Ratio 
Prediction  Results 
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see  to  /m  fork  raids 


Fig  53  Nozzle  3  Side  t.o  Axial  Force  Ratio 


Prediction  Results 


Prediction  Results 
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match  the  measured  data  trends.  On  the  other  hand,  the 
Green/McCullough  model  with  =  1  does  appear  to  provide 
some  correlation  with  the  measured  data.  This  model 
requires  an  empirically  derived  value  for  ,  which  was 
selected  by  iteration  to  be  the  value  that  gave  the  smallest 
average  error  from  measured  data  for  all  vectored  flow 
cases.  Frorr.  observation,  some  data  trends  seem  to  follow 
the  Green/McCullough  model,  except  their  values  were  too 
large.  For  these  cases  it  seemed  that  adding  an  empirical 
constant ,  .  to  the  Gr een/McCul 1 ough  equation  Cwith  =  ID 

might  provide  better  correlation  with  the  measured  data  than 
multiplying  by  the  coefficient  .  Again,  this  constant  was 
derived  by  iteration  to  produce  the  lowest  error  between 
predicted  and  measured  data.  Table  8  lists  the  resulting 
coefficients  CC^D  and  constants  CC^D  and  their  appropriate 
errors  for  each  nozzle.  Figures  55  through  58  graphically 
illustrate  the  correlation  between  the  measured  and  modeled 
data.  For  nozzles  1  and  2,  the  average  errors  for  using  CL, 
CIT.  3  and  10.9  percent!  were  significantly  less  than  the 
average  errors  for  using  .  For  nozzles  3  and  4,  the 

average  error  for  using  C26.6  and  28.1  percent}  were 

slightly  less  than  the  average  errors  using  C^- 

Shock  Loss  Model 

Although  the  Green/McCullough  equation  can  provide  an 
approximate  value  for  side  to  axial  force  ratio,  it  does  not 
provide  any  information  concerning  the  magnitude  of  the 
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Table  8. 


Average  Error  from  Actual  Results  for 


Gr een/McCul  1  ough  Side  to  Axial  Force  Prediction  Model 


Nozzl e 

Model 

Coef  f l ci ent 

Ci 

Error  using 

Model  with 

Coef  f i ci ent 

Model 

Constant 

C2 

Error  using 

Model  with 

Constant 

1 

4.  29 

29.  9 

0.  200 

17.3 

2  * 

3.  71 

18.3 

0.  218 

lO.  9 

3  * 

1 . 44 

26.  6 

0.  039 

28.  1 

4  * 

0.  69 

80.  7 

O.  029 

81 . 1 

»  data  for  non-vectored  cases  not  considered 


SEE  TO  AMAL  FORCE  RATO  SOE  TO  AXIAL  FORCE  RATOS 


NOZZLE  1  PREDICTION  MODELS 
MODFED  GREEN/McCULLOUGH  MODELS 


Fig  55  Nozzle  1  Side  to  Axial  Force  Ratio  Prediction 
Results  with  Modified  Green/McCullough  Model 

NOZZLE  2  PREDICTION  MODELS 
MODFIED  GREEN/McCULLOUGH  MODELS 


MASS  FLOW  RATO  (MFR) 


Fig  56  Nozzle  2  Side  to  Axial  Force  Ratio  Prediction 
Results  with  Modified  Green/McCullough  Model 
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SCE  TO  AXIAL  FORCE  RAIDS  SOE  TO  AXIAL  FORCE  RATD 


NOZZLE  3  PREDICTION  MODELS 
MODFED  GREEN/McCULLOUGH  MODELS 


MASS  FLOW  RATO  (MFR) 


Fig  57  Nozzle  3  Side  to  Axial  Force  Ratio  Prediction 
Results  with  Modified  Green /Me Cul 1 ough  Model 

NOZZLE  4  PREDICTION  MODELS 
MODFED  GREEN/McCULLOUGH  MODELS 


MASS  FLOW  RATO  (MFR) 


Fig  58  Nozzle  4  Side  to  Axial  Force  Ratio  Prediction 
Results  with  Modified  Green/McCullough  Model 


forces  that  the  nozzle  produces, 


Earlier  calculations  of 


thrust  efficiencies  in  section  IV  of  this  report  show  that 
ideal  thrust  equations  do  not  always  provide  accurate  values 
of  the  actual  nozzle  thrust.  In  order  to  be  able  to  predict 
the  nozzle  total  thrust,  the  author  has  developed  a  model 
which  calculates  a  loss  ratio  for  total  pressure  used  to 
determine  a  nozzle  exit  velocity  which  is  applied  to  the 
thrust  equation.  For  this  simple  model,  other  losses  such 
as  those  due  to  wall  friction  are  not  considered.  The 
following  paragraphs  describe  the  development  of  this  model. 

The  thrust,  equation  for  a  nozzle  can  be  written  as: 


F  =  Cm  -►  m  .“)  V  +  CP  -  PDA  C7D 

t  s  p  e  e  a  e 


where  F4  is  the  total  thrust,  m  and  m  are  the  secondary 
t  s  p 

and  primary  mass  flow  rates,  is  the  nozzle  flow  exit 

velocity,  P  and  P  are  the  nozzle  exit  and  ambient  static 
e  a 

pressures,  and  A  is  the  nozzle  exit  area  C15:49I>  In  order 

e 

to  account  for  the  thrust  loss  due  to  the  radial  divergence 
of  streamlines  in  the  nozzle,  the  value  of  must  be 
ad j  usted : 
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C  85 
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1  +  cos  a 
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where  X  Is  the  loss  coefficient  for  radial  divergence  and  a 
is  the  nozzle  half  angle  of  the  diverging  section,  which 


equals  20°  for  all  test  nozzles  for  this  report  C20:2335. 

For  this  report,  all  equations  except  were  measured  using 
relatively  simple  instrumentation.  Therefore,  the  goal  of 
the  author’s  model  is  to  predict  the  value  of  and  apply 
it  to  the  thrust  equation  in  order  to  predict  the  total 
thrust  of  the  nozzle. 

Standard  isentropic  equations  provide  an  ideal  exit 

velocity  for  nozzles  C19:2385.  However,  as  discussed  in 

section  IV  of  this  report,  the  nozzles  tested  suffered 

thrust  losses  due  to  internal  shocks.  If  a  total  pressure 

ratio  for  shock  loss  were  determined,  the  isentropic 

equation  for  exit  velocity  could  be  used  with  this  ratio  to 

predict  the  actual  values  of  V  .  This  ratio  will  be  defined 
f  e 

as : 
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SL 
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where  77^  is  the  total  pressure  ratio  for  shock  loss,  is 

the  nozzle  flow  total  pressure  before  shock  loss,  and  P  is 
the  nozzle  flow  total  pressure  after  shock  loss.  Using  this 
ratio  and  the  isentropic  exit  velocity  equation,  the 
equation  for  the  calculated  exit  velocity  including  shock 
loss  effects  becomes; 
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where  y  is  the  specific  heat  ratio,  R  is  the  universal  gas 
constant,  is  nozzle  flow  total  temperature,  P&  is  exit, 
static  pressure  and  P^  is  nozzle  supply  total  pressure. 

From  the  results  and  observations  of  sections  III  and  IV  of 
this  report,  it  is  reasonable  to  assume  that  W1  1 1  be 

related  to  the  pressure  losses  resulting  from  oblique 
shocks.  In  order  to  determine  the  total  pressure  loss 
across  an  oblique  shock,  the  Mach  number  upstream  of  the 
shock  and  the  shock  angle  must  be  known.  The  lsentropic 
Mach  number  at  the  secondary  flow  port  can  be  determined  by 
the  ratio  of  nozzle  area  at  the  port  to  the  nozzle  throat 
area.  Oblique  shock  angle  can  be  calculated  by  the 
equati on: 

p  h 

1  y  +  1  2  y  -  1 

ft  =  arc  sin  •  -  -  -  +  - 

«1  P1  J  ^ 


1  u.  J  « 

Cl 


where  ft  is  the  oblique  shock  angle  with  respect  to  the  flow 
turning  angle,  M^  is  the  Mach  number  before  the  shock,  y  is 
the  specific  heat  ratio,  P^  is  the  static  pressure  before 
the  shock,  and  P  is  the  static  pressure  after  the  shock 
C6:4633.  Figure  59  illustrates  how  the  value  of  P^  and  P^ 
can  be  taken  from  the  measured  values  at  static  ports  P9  and 
PI  2  of  the  test  nozzle. 
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=  static  pressure  upstream  of  oblique  shock  C=  PS) 

P^  =  static  pressure  downstream  of  oblique  shock  C=  P12D 


Fig  59.  Oblique  Shock  Formation  in  Nozzle  Diverging  Section 
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Once  the  oblique  shock  angle  has  been  determined,  . 
the  Mach  number  downstream  of  the  oblique  shock  can  be 
calculated  as: 


M 


2 


2  _ 

sin  ft 

tan  C  6 

-  m 

(  r  *  i  ' 

- 

r  ■ 1 1 1 

tan 

ft 

l  2  . 

l  2  J  J 

-0.  5 
Cl  3D 


where  ft  is  the  oblique  shock  angle.  6  is  the  flow  turning 
angle,  and  y  is  the  specific  heat  ratio  C19:36CQ.  For 
nozzle  flow  with  secondary  flow  applied  from  one  wall  ,  a 
value  of  6  =  40°  is  used  to  model  the  flow  turning  to  become 
parallel  to  the  opposite  nozzle  wall.  For  nozzle  flow  with 

O 

secondary  flows  applied  from  both  walls,  a  value  of  6  =  20 
is  used  to  model  the  flow  turning  to  become  parallel  to  the 
nozzle  axis.  From  the  values  of  M^  and  M^.  the  isentropic 
total  to  static  pressure  ratio  can  be  determined: 


C14D 


Using  these  ratios,  the  value  of  the  total  pressure  ratio 
due  to  shock  loss  is  simply: 


In  summary,  the  above  loss  model  requires  the  following 
i nputs : 

rrip  =  primary  mass  flow  rate 
m,^  =  secondary  mass  flow  rate 
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P  =  ambient  static  pressure 
a 

P  =  nozzle  exit  static  pressure 
e 

P^  =  nozzle  primary  flow  supply  pressure 

P^  =  nozzle  centerline  static  pressure  at  secondary 
f 1 ow  por  t 

P^  =  nozzle  centerline  static  pressure  downstream 
of  secondary  flow  port 

A  =  nozzle  exit  area 
e 

When  the  total  force  values  calculated  using  this  loss 
model  were  compared  to  measured  values,  the  model  values 
were  consistently  larger.  Therefore,  an  additional  loss 
coefficient  of  0.9  was  added  to  equation  8  for  the  case  of 
secondary  flow  being  applied  from  both  nozzle  walls.  This 
additional  total  pressure  loss  coefficient  reflects  the 
formation  of  Mach  reflections  instead  of  simple  oblique 
shocks  for  these  cases.  Table  9  lists  the  average  ratios  of 
the  total  force  calculated  using  the  shock  loss  model 
developed  above  to  the  measured  total  force.  From  this  one 
can  see  that  the  shock  loss  model  provides  an  excellent 
prediction  of  the  nozzle  total  force.  Figures  60  through  63 
graphically  illustrate  the  correlation  between  the  shock 
loss  model  total  forces  and  actual  measured  total  forces. 
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Table  9  Average  Ratio  of  Total  Force  Calculated  Using  the 


Shock  Loss  Model  to  Measured  Total  Force 


Ratio  of  Loss 

Model  to  Measured  Total  Force 

Nozzl e 

Secondary  Flow 
Applied  from 
Upper  Wall 

Secondary  Flow 
Appl i ed  f r  om 
Lower  Wall 

Secondary  Flow 
Applied  from 
Both  Walls 

1 

O.  952  ±  O.  058 

O. 952  ±  O. 057 

1  .  Ol  ±0.  008 

2 

O.  992  ±  O.  057 

0. 956  ±  O. 053 

1 . 01  tO.  033 

3 

O.  901  ±  0.  OlO 

O.  893  ±  O.  013 

O.  978  t  0.  018 

4 

O.  985  t  0.  017 

O.  995  i  0.  015 

O.  995  ±  0.  014 

lOl 


TOTAL  FORCE  (LBF)  TOTAL  FORCE  (L0F) 


NOZZLE  1  LOSS  MODEL  COMPARISON 
AVERAGE  NPR  =  42.2 
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Fig  60  Comparison  of  Loss  Model  and  Measured  Total 
Force  Values  for  Nozzle  1 


Fig  61  Comparison  of  Loss  Model  and  Measured  Total 
Force  Values  for  Nozzle  2 


1  02 


TOTAL  FORCE  (LflF)  |  TOTAL  FORCE  (LBF) 


Fag  62  Campai  l^on  of  Loss  Model  and  Measured  Total 
Force  Values  for  Nozzle  3 


Fig  63  Comparison  of  Loss  Model  and  Measured  Total 


Force  Values  for  Nozzle  4 
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VI .  Cone  I usi ons 


From  the  observations,  results,  am  analyses  presented 
in  this  report,  the  following  conclusions  can  be  made: 

1.  The  presence  of  secondar /  injection  ports  may  induce 
oblique  shocks  to  form  in  the  diverging  section  of  the 
nozzles  tested.  These  shocks  cause  the  flow  to  separate 
from  the  nozzle  wall.  The  interactions  of  these  shocks 
influence  the  flow  path  within  the  nozzle. 

2.  The  C.JTVC  nozzles  tested  showed  increased  axial  flow 
stability  with  increasing  NPR  for  primary  flows  with  no 
secondary  flow. 

3.  A  critical  secondary  flow  to  primary  flow  mass  flow 
ratio  CMFR3  exists  below  which  a  nozzle  wlII  not  vector.  A 
second  critical  MFR  exists  that  produces  maximum  vector 
performance  within  a  given  nozzle. 

4.  Flow  visualization  and  side  force  measurements  show 
that  all  nozzles  had  a  slight  bias  towards  flow  attachment 
to  the  lower  nozzl  e  wall,  therefore  producing  a  positive 
vector  angle  bias. 

5.  From  a  comparison  of  the  flow  performance  parameters 
measured  at  a  primary  flow  supply  pressure  of  2.923  bar 
absolute,  nozzle  1  has  the  best  vectoring  performance  of  all 
the  nozzles  tested. 

6.  Flow  vectoring  response  times  for  nozzle  1  at  its 
best  vector  range  MFR  vary  from  43.0  to  140.2  ms  C7.1  to 
£3.  3  Hzl>  when  measured  from  the  point  of  time  when  the 
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secondary  flow  solenoid  valve  opens  or  closes.  The  same 
response  times  vary  between  3. 0  and  98. 2  ms  CIO.  2  to  333.  3 
Hz3>  when  measured  from  the  point  in  time  when  the  secondary 
flow  reaches  its  supply  line  pressure. 

7.  The  Green/McCullough  model  for  side  force  to  axial 
f  or  ce  r  atio  provides  a  good  approximation  of  the  actual 
ratio  when  the  proper  multiplicative  coefficient  or  additive 
constant  is  used. 

8.  The  shock  loss  model  developed  in  this  report 
provides  an  excellent  prediction  of  the  actual  total  force 
produced  by  the  nozzles  tested. 
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VI I  . 


Recommendati ons 


For  anyone  desiring  to  continue  research  in  the  areas 
described  in  this  report,  the  following  recommendations  are 
pr  ovi  ded : 

1.  Replace  the  existing  plexiglass  walls  of  the  nozzle 
test  assembly  for  schlieren  photographs  with  opti cal -qual 1 ty 
glass.  This  would  allow  for  better  photographs  of  shock  and 
boundary  layer  structures  within  the  nozzle. 

2.  Study  nozzle  starting  process  at  lower  NPRs  to 
determine  causes  of  flow  bias,  especially  the  formation  of 
vortices.  High  speed  cameras  and  oil  film  or  china  clay 
flow  visualization  could  be  used  in  this  study. 

3.  Develop  permanently  sealed  nozzle  assemblies  rather 
than  interchangeable  nozzle  wails  and  plates.  This  will 
help  avoid  any  flow  bias  due  to  nozzle  leakage  of  assembly 
mi  sal l gnment 

4.  Study  the  effects  of  using  suction  in  the  secondary 
flow  por  i.  for  enhanced  flow  attachment  to  the  wall  opposite 
the  secondary  port  which  is  injecting  flow. 

5.  Study  the  effects  of  varying  the  angle  of  the 
secondary  flow  axis  to  where  it  points  upstream  into  the 
pr  l  mar  y  f  1  ow. 

0  Modify  the  existing  test  assemblies  to  move  the 
secondary  flow  solenoid  valves  as  close  as  possible  to  the 
secondary  injection  ports. 
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Appendi x  A:  Test  Model  and  Test  Facilities 
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Fig  64  Nozzle  Design  Details  and  Dimensions 


Fig  65  VKI  Wind  Tunnel  Apparatus 


109 


no 


Appendix  B: 


Sample  Calculations 

-  Thompson  Separation  Model 

-  Ideal  Thrust  and  Thrust  Efficiency 

-  Side  Momentum  Gain  CG..) 

M 

*  Flow  Momentum  Ratio  (R  } 

M 

-  Reynolds  Number 

-  Boundary  Layer  Thickness 

-  Broadweli  Force  Ratio  Model 

-  Gr een'McCul  1  ough  Force  Ratio  Model 

-  Shock  Loss  Total  Force  Model 
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Thompson  Separation  Model 


Thompson  developed  the  following  equation  for  predicting 
the  Mach  number  at  flow  separation  in  a 

convergent-divergent  nozzle  subjected  to  an  adverse  pressure 
gr  jidi  ent  C  20 .  4 ?  : 
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'Jsing  isenticpic  relations,  the  area  ratio  at  which  this 
MaO.  occurs  can  be  determined  C7:  80: 
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where-  A  is  the  nozzle  area  at  which  M  occurs  and  is  the 
nozzle  throat,  area  For  constant  channel  height  nozzle 
geometry,  the  following  equation  for  the  length  along  the 


nozzle  axis  at  which  separation  occurs  Cx  5  is: 


X 

- 

1 

Wt 

1  A  1 

'  Wt 

2  tan  o 

> 
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i _ 

SAMPLE  CALCULATION: 


Input  Values: 


Output  Values: 


£0 

0.050  C nor  mall  zed  va) ue5 

3.  09 

4.  61 


Cl  85 
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C. 248  C nor  mall  zed  value) 


x 

s 

Ideal  Thrust  and  Thrust  Efficienc, 

For  ideal  thrust  calculations,  the  flow  is  assumed  to 
1  sentr  opi  cal  1  y  expand  to  where  the  exit  static  pressure 
equals  the  ambient  static  pressure  08:26).  This  reduces 
equation  7  of  this  report  to: 


F 


t , i deal 


i  *)  V 

P  J  e 


0  9) 


where  V^  is  calculated  using  equation  11  of  this  repo;  t  with 

'n._,  =  1  C  no  shock  loss)  and  P  =  P  .  Thrust  efficiency  is 
SL  e  a 

simply  the  ratio  of  the  measured  total  thrust  to  the  ideal 
total  thrust: 


THRUST  EFFICIENCY 


^t , i deal 


SAMPLE  CALCULATION: 
Input  Values.  F 

nS 


Output  Values: 


m 

P 


P 

e 


r 


R 


V 

e 


=  0.  603  lbf 

=  0.000971  lbm/s 

=  0.  01 1 1  62  1  bm/s 

=  530  °R 

=  42.  4  psi a 

=  P  =113  psi  a 
a 

1  .  4 

1716  ft2/s2  °R 
=  2026  ft/s 

,  =  0,  764  lbf 

i  deal 


C  20) 
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THRUST  EFFICIENCY 


0.  789 


Side  Momentum  Gain 

The  side  momentum  gain  CGO  is  defined  as  the  ratio  of 

M 

measured  side  force  CF^)  to  the  momentum  of  the  secondary 
jet  flow,  and  is  given  by: 


G 


M 


C21)> 


where  m  is  the  secondary  mass  flow  rate  and  V  is  the 
s  s 

secondary  flow  velocity  V  is  calculated  using  equation  11 
of  this  report  with  7),,^  =  1  C  no  shock  loss)  ,  with  P^  equal 
to  the  secondary  supply  pressure,  and  with  equal  to  the 
static  pressure  measured  at  the  exit  cf  the  secondary  flow 
por  t . 


SAMPLE  CALCULATION: 


Input  Values; 

F 

s 

=  0.  19b  1  bf 

m 

s 

=  0.000971  lbm/s 

o 

T  = 

530  R 

Pt 

=  17.2  psi a 

P 

e 

=  4 .  22  psi a 

r  = 

1  .  4 

R  = 

P  2  o 

1716  ft  w/s  R 

Output  Values: 

V 

s 

=  1451  ft/s 

°M 

ii 

A 

A 

'J 
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Flow  Momentum  Ratio 


The 

seconder 

secondar 


Mow  momentum  ratio  1  R.O  is  defined  as  the 

M 

y  flow  momentum  to  primary  l  low  momentum  at 
y  flow  por  t ,  and  is  given  by: 

m  V 


R 


M 


m  V 


P  P 


ratio  of 
the 


C225 


wh^r  e  m__  and  are  th  '•  same  quantities  as  in  equation  21  . 

m  is  primary  mass  flow  rate,  and  V  is  the  flow  velocity  at 
P  P 

the  secondary  port.  is  given  by  the  lsentropic 

equa  1 1 ons : 


C23D 


C241 


whe r e  M1  .  T.  are  the  Mach  number,  static  temperature,  and 

total  temperature  at  the  secondary  port  C 7 : 551 .  M  is 

oo 

calculated  using  equation  17  of  this  report. 


sample  calculation: 


1 nput  Val  ues : 


2.  84 


m  =  0.  000971  lbm/s 


m  =  0.  01 1 1  62  1  bra/s 
P 

T  =  530  °R 


P 


t 


6#cond»r y 


17.2  psi  a 
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4 . 22  psi  a 


Output  Values: 


e , sec  ondar  y 
/  -1.4 

P  =  1716  ft2/s2  °R 

V  =  1451  ft/s 
s 

M  =  2  58 

-  1907  ft 7s 

Pw  -  0.  662 

M 


lie 


Reynolds  Number 

The  length-based  and  hydraulic  diameter -based  Reynolds 
numbers  are  defined  as: 


Re 


L 


u  L 
00 


V 


css: 


Re. 


u  D. 
on  h 


265 


wher  ■-  L  is  the  length  alcna  a  nozzle  wall  ,  u  is  f  low 
velocity,  i'  is  flow  kinematic  viscosity,  and  is  hydraulic 
diameter  defined  by 


D 


h 


4  A 
P 


C  275 


where  A  is  the  area  oi  the  flow  channel  and  P  is  the 
perimeter  of  the  flow  channel.  If  the  value  of  Re  and  Re 

h 

s 

are  greater  than  3  6  x  10  ,  the  flow  is  considered  to  be 
t  ur  bu  1  ent  fl  6:  415. 


SAMPLE  CALCULATION: 

I nput  Va 1 ues :  L  -  0 .  239  in  =  O.  Ol  992  f  t 

A  =  0.0619  in2  =  4.30  X  10"4  ft2 
P  =  O  996  in  =  0.0830  ft 
u  =  1907  ft/s 

CD 

v  =  3.  48  x  lo’5  f  t2/s 

Output  Values:  L>h  =  0.249  in  =  0.021  ft 
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TURBULENT  FLOW 


Re 


L 


1.09  x  10 
=  1.15  x  10h 


TURBULENT  FLOW 


Boundary  Layer  Thickness 

Equation  1  of  this  report  is  used  to  calculate  the 
boundary  layer  thickness  at  the  secondary  port  location. 
Equations  2  and  3  are  then  used  to  determine  displacement 
thickness  and  momentum  thickness. 


SAMPLE  CALCULATION 


I  nput 

Val ues : 

x  - 

0.  225  in  = 

O.  01992  ft 

u 

ao 

-  1907  ft/s 

v  - 

3  48  x  10~5 

2 

ft  /s 

Output 

Val  ues : 

\ 

=  0.  00548  in 

61 

=  0 . 00059  l  n 

62 

=  0.  00053  l  n 

0. 00053  l n 


E<r u.idwel J  Force  Ratio  Model 

Equation  4  of  this  report  is  used  to  calculate  the 

Broadwel  1  model  values  for  side  force  to  axial  force  ratio. 

is  calculated  using  equation  11  of  this  report  with 

n  ,,  =1  t  no  shock  lose)  and  P  =  P  .  V  is  calculated  usina 

t>L  e  a  s 

the  same  method  as  described  for  G  calculations  Cequation 

M 

£11  of  tins  :  eporO. 


SAMP!  .F  CALC'JLATI  ON  . 


i  fiput  Va  I  ue>_; 


}■  -  1-4 
c*  -  0  17 


M  =  2.  56 

at. 

V  =  1  907  f  t  's 

UJ 

P.  .  =  17  .  2  psi  a 

t . secondar  y 


e , secondar  v 


=  4 . 22  ps l a 


P,  =  42.  4  psi a 

t.  ,  pr  l  mar  y 

P  -  P  =  1  .  1 3  psi  a 

e ,  pr  i  mar  y  a 

P  =  1716  ft2/s2  °R 


T.  =  T,  -  530  R 
1 1  too 

m  -  O.  000971  I  bm/s 

l 

m  =  0.011162  1  bm/s 
P 


a 


90 
n  =  70 ' 


Output  Values: 


V  -  2026  ft/s 
ex 


1451  ft/s 

0.  241 


1  1  9 


;CuX  lough 


Equation  S  of  this  report  Is  used  to  calculate  the 
Green/"'McCul  lough  model  values  for  side  force  to  axial  force 
ratio.  I  is  calculated  using  equation  6  of  this  report. 
Two  values  of  the  equation  are  given;  one  value  for  when  the 
coefficient  is  multiplied  to  equation  5  and  one  value  for 
when  the  constant  Cg  is  added  to  equation  9  Cwith  *  ID. 


SAMPLE  CALCULATION: 


I nput  Val ues 


Output  Val ues : 


Cx  *  4.  29 
C2  =  0.  200 
wA  «  0.  000971  lbf/s 
w  «  0.011162  lbf/s 


P 

F  *  0.  603  lbf 


P 


t » secondary 


1.187  bar  abs 


**e .  secondar  y 
I  *  40.7  s 


O.  291  bar  abs 


*=  0.  175 


*  0,  241 
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Shock 


Loss  Total  Force  Model 
Squat-ions  7  through  15  of  this  report  are  used  to 
calculate  the  total  force  value  using  the  shock  loss 
modeling  developed  by  the  author.  The  order  of  calculation 
is  as  follows: 


1 . 

Use 

equation 

12 

to 

cal culate 

ft 

2. 

Use 

equation 

13 

to 

calculate 

M2 

3. 

Use 

equation 

14 

to 

calculate 

Pt/  P 

4. 

Use 

equation 

15 

to 

calculate 

^SL 

5. 

Use 

equation 

11 

to 

calculate 

V 

e 

6. 

Use 

equation 

9  to  calculate  \ 

7. 

Use 

equations  7 

and  8  to  calculate  F 

SAMPLE  CALCULATION: 

Input  Values:  =  2.58 

P2  *  5.  8©  psia 
P1  «  2.  21  psia 
r  =  1.4 
<5  »  40* 
o  =  20* 

m  *  0.000971  lbm/s 
s 

m  “  0.  011162  lbm/s 
P 

P#  «  3. 39  psia 
P^  ■  1.13  psia 

A  =  0.08397  in2 

e 

Output  Values:  ft  *  37.9* 

Mg  =  0.617 
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O.  184 


V  =  1162  ft/s 

© 

X  =0.97 
F  =  0.  609  lbf 
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The  nozzle  mass  flow  rates  were  measured  using  a  model 

01.1000  ROTAMETER  which  consists  of  a  wedge-shaped  float  in 

a  tapered  glass  tube  with  an  inlet  fitting  on  the  bottom  and 

an  exit  fitting  on  the  top.  The  principle  of  operation  of 

the  device  is  that  the  equilibrium  between  the  float’s 

weight,  buoyancy,  and  drag  will  place  it  at  a  certain  height 

in  the  tube  for  a  given  mass  flow  rate  and  pressure.  The 

values  of  float  height  and  fluid  pressure  are  recorded  and 

the  mass  flow  rate  is  calculated  using  formulas  contained  in 

ROTAMETER  company  literature.  A  metal  float  with  a  mass  of 

3 

59.3  g,  density  of  7.94  g/cm  ,  and  maximum  diameter  of  81.0 
mm  was  used  to  measure  primary  flow.  A  plastic  float  with  a 

3 

mass  of  15.28  g,  density  of  2.140  g/cm  .  and  maximum 
diameter  of  21.8  mm  was  used  to  measure  secondary  flow. 

Three  measurements  were  recorded  and  averaged  for  each  flow 
condition.  Fig  83  summarizes  the  mass  flow  rate  measurement 
results  and  compares  it  to  isentropic  mass  flow  rates 
calculated  assuming  choked  flow  at  the  nozzle  throat. 
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MASS  FlOW  flJN/S) 


MASS  FLOW  VALUES 


ROTAMETER  VS  SENTROPC  VALUES 


MASS  FLOW 
CALCULATIONS 

DIAL  TRANS 
PRESS  PRESS 
(BAR  ABS) (BAR  ABS) 

DIM.  ROTA  I SEN 

ERROR  NASS  FIONNASS  FLOW 
U>  (LBN/S)  (LBM/S) 

I  SEN 
ERROR 

(X) 

ROTA  HT 

(CM) 

PRIMARY  LINE 

1.20 

1.182 

-1.50  0.002366  0.005916 

150.01 

1.92 

1.50 

1.489 

-0.73  0.004346  0.007452 

71.49 

3.32 

2.00 

1.982 

-0.90  0.006964  0.00992 

42.46 

4.32 

2.50 

2.446 

-2.16  0.008973  0.012242 

36.43 

4.97 

3.00 

2.923 

-2.57  0.011162  0.01463 

31.07 

5.55 

3.50 

3.439 

-1.74  0.01284  0.017212 

34.05 

6.07 

+1  SECONDARY 

1.20 

1.187 

-1.08  0.000971  0.001903 

96.03 

1.60 

1.36 

1.284 

-5.59  0.001292  0.002059 

59.40 

2.05 

1.51 

1.447 

-4.17  0.001645  0.00232 

41.02 

2.55 

1.81 

1.685 

-6.91  0.002152  0.002702 

25.55 

3.12 

2.11 

1.976 

-4.35  0.00268  0.003169 

18.24 

3.55 

2.41 

2.236 

-7.22  0.003253  0.003585 

10.21 

3.80 

-1  SECONDARY 

1.20 

1.183 

-1.42  0.000969  0.001897 

95.70 

1.60 

1.36 

1.262 

-7.21  0.001281  0.002024 

58.03 

2.05 

1.51 

1.436 

-4.90  0.001639  0.002303 

40.48 

2.55 

1.81 

1.675 

-7.46  0.002146  0.002686 

25.18 

3.10 

2.11 

1.955 

-7.35  0.002666  0.003135 

17.61 

3.52 

2.41 

2.209 

-8.34  0.003234  0.003542 

9.54 

3.80 

Fig  07  Nozzle  Mass  Flow  Measurements 
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Appendix  D 


The  figures  on  the  following  pages  are  the  individual 
static  pressure  profiles  used  to  generate  the  summary  plots 


of  centerline  pressures  found  in  this  report.  Pressure 
values  shown  in  the  nozzle  outlines  are  given  in  bar 
absolute.  The  following  is  an  index  for  determining  which 
individual  figures  match  which  summary  figures: 


Summary  Figure  Nozzle  Individual  Figure 


Fig 

23 

0 

1 

3 

4 

5 

08 

09 

70 

71 

72 

Fig 

44 

1 

73 

2 

70 

3 

79 

4 

82 

Fig 

45 

1 

2 

3 

4 

74 

77 

80 

83 

Fig 

40 

1 

75 

2 

78 

3 

81 

4 

84 
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TEST  NUMBER  NPR  =  4-2.16  NPR  (CALC)  =40.80 


890131.21  +1  SI  = 

NOZZLE  1  -1  SI  = 

4-1 1  SI 


PORT  40  =  0.075 


iJ.227 

0.24Z- 

0.221 

0.221 

0.235 

0.272 

0.188 

0.215 

0.252 

0.281 

o.: 

0.229 

0.235 

0.260 

0.297 

9^ 

TOTAL  FORCE  =0.236 

AXIAL  FORCE  =0.234  VECTOR  ANGLE  =  6.2 
SDE  FORCE  =  0.025 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890131.21 


UPPER  WALL. 
UPPER  CVTY 
CENTERLtC 
LOWER  CVTY 
LOWER  WALL 
GEOCTRY 


ai  Qj  as  Q7  ae 

NORMALIZED  06TANCE  FROM  NOZZLE  TVROAT 


Fig  ©Q  Static  Pr»ssur»  Profiles  for  Test  Number  800131.21 


129 


TEST  NUMBER  NPR  =  39.30 

890201.39  +1  SI  =  0 

NOZZLE  3  -1  SI  =  0 


NPR  (CALC)  =37.18 
PORT  40  =  0.081 


+1i  SI 


3224 


1228  0.2^ 


0.220  0-224  0-245  0.279  0> 

0.1 83  0.224  0.259  0.282  0. 

0-224  0.226  0.251  0-290 


-I1  SI 


TOTAL  FORCE  =0.238  I 

AXIAL  FORCE  =0.237  VECTOR  ANGLE  =  4.7 
SIDE  FORCE  =  0.019 


Fig  70  Static  Pratsur#  Profile®  for  Tost.  Number  80C201 . 30 


130 


:  (BW  ABSOLUTE) 


TEST  NUMBER  NPR  =  40.16  NPR  (CALC)  =37.26 

890207.06  +1  SI  =  0  PORT  40  =  0.079 

NOZZLE  4  -1  SI  =  0 


-1i  SI 


5034  0.050 


0.052 


0.024  0-032  0-042  0.050  O.C 

1  ^  ^050  °-040  °-036  0-044  0.052  O.C 

0.034  0.044  0  048  O.C 

^  Sl  Q.0fi3  0-050 

TOTAL  FORCE  =0.238 

AXIAL  FORCE  =0.238  VECTOR  ANGLE  =  3.1 
SIDE  FORCE  =  0.013 


Fig  71  Static  Pressure  Profiles  for  Test  Number  800207.00 
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PtESSlK:  (BAR  IWSCIUTE) 


-0.2 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890126.32 


0.1  03  '  as  07  OB 

NORMALIZED  DGTANCC  FROM  NOZZLE  TVWOAT 


Fig  72  Static  Praissur*  Profile  Tor  T*sL  Number  800126.  32 
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Fig  73  Static  Pressure  Profiles  for  Test  Number  800131.25 
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NPR  (CALC)  =40.26 
PORT  40  =  0.077 


Fig  74  Static  Pressure  Profiles  for  Test  Number  800131.28 
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PHESSLFE  (BAR  ABSOLUTE) 


TEST  NUMBER  NPR  =  40.38  NPR  (CALC)  =39.47 

890131.27  +1  SI  =  1.2  PORT  40  =  0.079 

NOZZLE  1  -1  SI  =  1.2 


+1i  SI 


5803  0.51 


0.483  0,503  0  511  0.509  o!1 

1 6A34%43(?.392  0,481  0,507  0,51 2  0,503  a; 

5;ZS;S^5g4^^4B7  0.517  0.524  °-5110, 

TOTAL  FORCE  =0.293 

AXIAL  FORCE  =  0,292  VECTOR  ANGLE  =  4.1 
SDE  FORCE  =  0.021 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890131.27 


UPPER  WALL 


UPPER  CVTY 


CENTERLf>C 


LOWER  CVTY 


LOWER  WALL 


GEOMETRY 


NORMALIZED  OSTANCE  PROM  NOZZLE  TVROAT 


Fig  75  Static  Pressure  Profiles  for  Test  Number  800131.27 
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PRESSURE  (BAR  AflSCH/TI) 


TEST  NUMBER  NPR  =  41.60  NPR  (CALC)  =40.56 

890202.48  SI  =  1.81  PORT  40  =  0.074 

NOZZLE  2  -1  SI  =  O 


4-1  V  SI 


0.442 

"0452 

0.446 

O.AS^ 

0.452 

0.464  0. 

376 

0.435 

0.450 

0.458 

o 

■4 

o 

1407 

0.448 

0.476 

0.491 

0.491  n, 

~1  Sl  ossa 

TOTAL  FORCE  =0.290 

AXIAL  FORCE  =0.280  VECTOR  ANGLE  =  16.6 
SDE  FORCE  =  0.080 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890202.48 


UPPER  WALL 
UPPER  CVTY 
CENTERLMi 
LOWER  CVTY 
LOWER  WALL 
GEOMTTRY 


02  QA  08  OB 

0.1  OJ  OB  07  as 

NORMALIZED  OSTANCE  FROM  NOZZLE  TVROAT 


Pig  7B  Sifttic  PrMiur*  Profil**  for  TmC  Number  800303.  48 
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TEST  NUMBER  NPR  =  40.31  NPR  (CALC)  =39.32 

890202.49  4-1  SI  =  0  PORT  40  =  0.077 

NOZZLE  2  -1  SI  =  1.81 


4-1  \  SI 


*472 


.538  0.5S 


0.448  0,479  0,493  0.493  0^37 

1^4af7i19&4^392  0,438  0,458  0466  °-470  0.3  ■ 47 

0.450  0.460  0488^ 

TOTAL  FORCE  =0.290 

AXIAL  FORCE  =  0.290  VECTOR  ANGLE  =-1 3.2 
SDE  FORCE  =  -0.070 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890202.49 


UPPER  WAUL 


UPPER  CVTY 


CENTERLNE 


LOWER  CVTY 


LOWER  WALL 


GEOC7RY 


02  O*  06  08 

1  03  05  07  08 

NORMALIZED  distance  from  nozzle  tvroat 


Fig  77  Static  Pressure  Profiles  for  Test  Number  860202. 46 


137 


T 


TEST  NUMBER  NPR  =  39.38  NPR  (CALC)  =39.02 

890202.50  +1  SI  =  1.81  PORT  40  =  0.078 

NOZZLE  2  -1  SI  =  1.81 


+1  v  SI 


\A2^8. 


’'’0624 

0.630* 

r<70 

[  0.577 

0.612 

0.616 

0.599  0. 

f.505  °-571 

0.603 

0.608 

0.589  0.- 

5  0.571 

mOI 

0.601 

0.608 

0.595  o 
y 

_1  Si  n  n 

TOTAL  FORCE  =0.340 

AXIAL  FORCE  =0.340  VECTOR  ANGLE  =~0.8 
SDE  FORCE  =-4.760 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890202.50 


— 

1 1  i 

n 

. 

■  V*VI 

_ - 

* 

— 

_ 3 

_ 

_ 

—— — 

tz3 

UPPER  WALL 


UPPER  CVTY 


LOWER  CVTY 


LOWER  WALL 


CEOMETRY 


Fig  78  Static  PrtEsur*  Profiles  for  T#st  Number  800202.50 
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(BAR  ABSOLUTE) 


0.23] 


NPR  = 

39.80 

NPR  (CALC) ; 

=  35.72 

-hi  SI  = 

2.41 

PORT 

40  = 

0.079 

-1  SI  = 

0 

i 

SI 

'755.37 

051^ 

0.520 

0.537 

0.514 

0.539  OT 

M 

19S,46§.481 

0.529 

0.555 

0.555 

0.549  0.3 

0.537 

0.570 

0.594 

0.580  H2 

m 

_1  Sl  ^ - -  0-S3f 

TOTAL  FORCE  =0.320  ! 

AXIAL  FORCE  =0.320  VECTOR  ANGLE  =11.1 
SIDE  FORCE  =  0.060 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890202.09 


UPPER  WAU. 
UPPER  CVTY 
CENTERLfC 
LOWER  CVTY 
LOWER  WALL 
GEOMETRY 


02  04  ae  ob 

Ol  03  as  07  Ofl 

NORMALIZED  DISTANCE  FROM  NOZZLE  TVROAT 


Fig  79  Static  Pressure  Profiles  for  Test  Number  890202.09 
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TEST  NUMBER  NPR  =  38.02 

890202.10  -4-1  SI  =  0 

NOZZLE  3  -1  SI  =  2.41 


NPR  (CALC)  =36.32 
PORT  40  =  0.083 


+1l  SI 


U oS'4^ 


^-0^86 

lK627 

0.643- 
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Fig  80  Static  Pratsur*  Profilat  for  Tact  Number  800202.10 
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PRES9JRE  (BAR  absolute) 


NOZZLE  PRESSURE  TEST  (CATON) 
TEST  NUMBER  890202.11 


Fig  81  Static  Pressure  Profiles  for  Test  Number  890202.11 
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Pig  8 2  Static  Prassur*  Profilaia  Tor  Tact  Numbar  800207.20 
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PRESSLRE  (BAR  ABSOLUTE) 
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Fig  84  Static  Pressure  Profiles  for  Test  Number  800207.22 


144 


B1  blioaraphv 


1.  Blondeaux,  P.  and  B.  de  Bernardinis  “On  the  Formation  of 
Vortex  Pairs  Near  Orifices,”  Journal  of  Fluid  Mechanics. 

1983. 

2.  Brins  made,  A.  "Thrust  Vector  Control,"  Solid  Rocket 
Technology,  edited  by  M.  Shorr  and  A.  Zaehringer.  New  York: 
Wiley  and  Sons,  Inc.  .  1967. 

3.  Broadwell  .  James  E.  “Analysis  of  the  Fluid  Mechanics  of 
Secondary  Injection  for  Thrust  Vector  Control," 

AIAA  Journal.  1:  1067-1075.  May  1963. 

4.  Carter,  J.M.  "Thrust-Axis  Control  In  Sol  id -Propel  1  ant 
Rocket  Motors,"  Rocket  Propulsion  Technology  Volume  1, 
Proceeding  of  the  First  Rocket  Propulsion  Symposium.  1-15. 
New  York.  Plenum  Press,  1961. 

5.  Cates.  James  W.  Two-Dimensional  Confined  Jet  Thrust 
Vector  Control  with  Flow  Visualization  and  Variable  Flow 
Geometry.  MS  Thesis  AFIT/GAE/AA/85D-2.  School  of 
Engineering,  Air  Force  Institute  of  Technology  CALD, 

Wright -Patter son  AFB  OH ,  December  1985. 

6.  Chapman,  Alan  J.  and  William  F.  Walker  Introductory 
Gas  Dynamics.  New  York:  Holt,  Rinehart,  and  Winston,  Inc.  , 
1971. 

7.  Clemow.  T.  Short-Range  Guided  Weapons.  London:  Temple 
Press,  1961. 

B.  Fitzgerald,  R.  E.  and  R.  F.  Karnpe  "Boundary  Layer  TVC  for 
Missile  Applications."  AI  AA/SAC^ASME  19th  Joint  Propulsion 
Conference  AI AA-83-1153.  American  Institute  of  Aeronautics 
and  Astronautics,  New  York,  1983. 

9.  Green,  C.  J.  and  Foy  McCullough  Jr.  "Liquid  Injection 
Thrust  Vector  Control,"  AIAA  Journal .  1 :  573-578,  March 
1963. 

10.  Kayes.  W.  M.  and  M.  E.  Crawford.  Convective  Heat  and 
Mass  Transfer.  New  York:  McGraw-Hill  Book  Company,  1980. 

11.  Lambert,  P.  G.  and  M.  E.  Franke.  "Characteristics  of 
Axisymmetric  Confined  Jet  Thrust  Vector  Control,"  AI AA/ASME/ 
SAE^ASEE  24th  Joint  Propulsion  Conference.  AI AA- 88 -3237 , 
Washington,  D.  C.  ,  1988. 

12.  Porzio,  Anthony  J.  "Expert mental  Study  of  a  Confined 
Jet  Thrust  Vector  Control  Nozzle,"  AIAA  24th  Aerospace 
Sciences  Meeting.  AI AA-86-0448 ,  Reno,  Nevada  1986. 


145 


13.  Schlichting,  Hsrmann.  Boundary-Layer  Theory.  New 
York:  irfcGraw-Hi.il  Book  Company,  1979. 


14.  Sharplo,  A.  H.  The  Dynamics  and  Thermodynamics  of 
Compressible  Fluid  Flow.  Volume  I.  New  York:  The  Ronald 
Press,  1953. 

15.  Sutton.  George  P.  Rocket  Propulsion  Elements.  Fifth 

•Edl tl on.  New  York:  John  Wiley  and  Sons.  1986. 

16.  Talda.  Timothy  A.  Experimental  Inyestiaation  of 
Two-Dimensional  Confined  Thrust  Vector  Control.  MS  Thesi s 
AF1  T''GAE^AAx'870-22.  School  of  Engineering,  Air  Force 
Institute  of  Technology  CALD,  Wri ght-Patter son  AFB  OH. 

Dec  ember  1 987 . 

17  Thompson.  R.  V.  Theoretical  Determination  of  the  Point 
of  Separation  in  a  Compressible.  Turbulent.  Supersonic 
Boundary  Layer  Subjected  to  an  Adverse  Pressure  Gradient. 
ASKL  Paper  67-WA/FE-38,  American  Society  of  Mechanical 
Engineers.  Maw  York,  1967. 

18.  von  Kantian  Institute  for  Fluid  Dyna~‘cs.  Facilities 
and  Instrumentation  1986.  Wyckmans  and  -mpany,  Brussels, 
Belgium.  198G. 

19.  Zucrow,  Maurice  J.  and  Joe  D.  Hoffman  Gas  Dynamics. 
Volume  I.  New  York:  John  Wiley  and  Sons,  1976. 

20.  Zukoskl .  Edward  E.  and  Frank  W.  Spaid.  "Secondary 
Injection  of  Gases  into  a  Supersonic  Flow,  "  AIAA  Journal  .  2: 
1089-1696.  October  1964. 


146 


Vita 


Captain  Jeffrey  L.  Caton  was  born  on  9  September  1960  in 
Arlington.  Virginia.  He  graduated  from  high  school  in 
Nokesville.  Virginia  in  1978  and  attended  the  University  of 
Virginia,  from  which  he  received  the  degree  of  Bachelor  of 
Science  in  Chemical  Engineering  in  May  1982.  Upon 
graduation,  he  received  a  commission  in  the  USAF  through  the 
ROTC  program.  He  reported  to  Wright -Patter son  AFB,  Ohio  in 
September  1982  where  he  completed  USAF  System  Acquisition 
Training  and  Reliability  and  Maintainability  training  in 
April  1983.  He  then  served  as  Lead  Product  Assurance 
Engineer  for  the  F/FB-111  Avionics  Modernization  Program. 
Deputy  for  Strategic  Systems,  Aeronautical  Systems  Division 
until  September  1084  when  he  became  Mission  Planning  Systems 
Integration  Engineer  for  the  Advanced  Cruise  Missile 
Program.  Deputy  for  Strategic  Systems,  Aeronautical  Systems 
Division.  Following  this  assignment,  he  entered  the  School 
of  Engineering,  Air  Force  Institute  of  Technology,  in  June 
1987. 


147 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


URITY  classification  authority 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  markings 


Form  Approved 
OMB  N o.  0704-0198 


I  I 


2b.  DECLASSIFICATION/ DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBFR(S) 

AFIT/GAE/AA/8  9M- 1 


6b  OFFICE  SYMBOL 
(If  applicable) 

AFIT/I24G 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

School  of  Engineering 


6c  ADDRESS  (City,  State,  and  ZIP  Cod*) 

Air  Force  Institute  of  Technology  (AU) 
Wright-Pattersan  AFB,  Ohio  45433-6583 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

von  Harman  Institute 


8b  OFFICE  SYMBOL 
Of  applicable) 


3.  DISTRIBUTION /AVAILABILITY  OF  REPOR 

Approved  for  public  release; 
Distribution  unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 

von  Harman  Institute  for  Fluid  Dynamics 
Chaus  see  de  Waterloo,  72 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

TASK 

NO 

NO 

WORK  UNIT 
ACCESSION  NO 


1 1 .  TITLE  (Include  Security  Claudication) 

TVO-DIMENSICNAL  CONFINED  JET  THRUST  VECTOR  CONTROL:  OPERATING  MECHANISMS  AND  PERFORMANCE 


Li 


12.  PERSONAL  AUTUORCj 

ffrey  L.  Caton,  Captain,  USAF 


PE  OF  REPORT  13b  TIME  COVERED 

Thesis  from _ to 


i6.  supplementary  notation 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

March,  1989  165 


17.  COSATI  COOES _ I  18.  SUBJECT  TERMS  ( Continue  on  reverie  if  neceuary  and  identify  by  block  number) 

FIELD  I  GROUP  I  SUB-GROUP  « _ *  r- _ j  —  _ 


Thrust  Control , Secondary  Injection,  Two-Dimensional, 
Thrust  Vector  Control  Systems, Flow  Visualization 


19.  ABSTRACT  ( Continue  on  reverie  if  neceuary  and  identify  by  blocl.  number) 

Dr.  Milton  E.  Franke 


ZC.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

unclassified/unlimitec  □  same  as  rpt  □  otic  users  Unclassified 


NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

Milton  E.  Franke,  Professor  513-255-2562  AFIT/ENY 


DDForm  1473,  JUN  86  Previous  editions  are  obsolete  _ security  classification  of  rn,$ 


SECURITY  CLASSIFICATION  OF  [H.S  PAGE 


Unclassif  red 


UNCLASSIFIED 

An  experimental  investigation  of  two-dimensional 
confined  jet  thrust  vector  control  nozzles  was  performed. 
Thrust  vector  control  was  accomplished  by  using  secondary 
flow  injection  in  the  diverging  section  of  the  nozzle. 
Schlieren  photographs  and  video  tapes  were  used  to  study 
flow  separation  and  internal  shock  structures.  Nozzle 
performance  parameters  were  determined  for  nozzle  flow  with 
and  without  secondary  flows.  These  parameters  included 
nozzles  forces,  vector  angles,  thrust' efficiencies ,  and  flow 
switching  response  times.  Vector  angles  as  great  as  18 
degrees  with  thrust  efficiencies  of  0.79  were  measured. 
Several  confined  jet  nozzles  with  variations  in  secondary 
flow  port  design  were  tested  and  results  were  compared  to 
each  other.  Converging-diverging  nozzles  of  similar  design 
to  the  confined  jet  nozzles  were  also  tested  and  results 
were  compared  to  the  confined  jet  nozzle  results. 

Existing  prediction  models  for  nozzle  side  to  axial 
force  ratio  were  evaluated.  A  model  for  nozzle  total  forces 
based  on  shock  losses  that  predicted  values  very  close  to 
actual  results  was  developed.  j-  ■  t.  <Jn j-, -vnn  ,*  *  /  i-  r'> 

»  ,  •  C*.  .  *  "frr  (  t  f  I  £  .  C  f-  >  -j .  -  ,  ,  % ,  . 

\  '  ^  ■ 

•'  '  ?  r  •  k'  •  •••  •<-,  !  :  <  -•  ‘ 


UNCLASSIFIED 


